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Partonic Matter Created in Au+Au Collisions 

at RHIC 

by 

Paul Richard Sorensen 

Doctor of Philosophy in Physics 
University of California, Los Angeles, 2003 
Professor Huan Z. Huang, Chair 

Measurements of identified particles over a broad transverse momentum px range 
may provide particularly strong evidence for the existence of a thermalized par- 
tonic state in heavy-ion collisions {i.e. a quark-gluon plasma). Of particular 
interest are the centrality dependence and the azimuthal anisotropy in the yield 
of baryons and mesons at intermediate px- The first measurements of v<i — an 
event-by-event azimuthal anisotropy parameter — and the nuclear modification 
factor R C p for mid-rapidity Kg and A + A production in Au+Au collisions at 
ultra-relativistic energy are presented. The K° S) A, and A candidates are selected 
based on characteristics of their decays in the STAR Time Projection Chamber 
(TPC). A statistical treatment is used to extract ^(pt) and Rcp(pt) from their 
invariant mass distributions. These measurements establish the particle type 
dependence of t>2 and Rep in the kinematic region 0.4 < px < 6.0 and \y\ < 1.0. 

In the low p T region (p T < 1.0 GeV/c) the v 2 values for different particles 
are increasing with p T and follow a mass dependence similar to that expected 



xvii 



from hydrodynamical models of Au+Au collisions — where, at a given p?, the 
particle with the larger mass will have a smaller v 2 . At higher pt however, v 2 
of the heavier A hyperon continues to increase while v 2 of the lighter K s meson 
saturates at v 2 ~ 0.13 for 2.0 < pr < 5.0 GeV/c. At intermediate pr the v 2 of 
K s and A + A are shown to follow a number-of-constituent-quark scaling with 

v% aon (p T /2) v l 2 ambda (p T /3) 
2 ~ 3 • 

The binary collision scaled centrality ratio R C p shows that A + A production 
at intermediate pt increases more rapidly with system size than kaon production: 
This is consistent with a scenario where multi-parton dynamics play an important 
role in particle production. At p T w 5.5 GeV/c A + A, K° S) and charged hadron 
production are all suppressed by a similar amount: a factor of three below ex- 
pectations from binary nucleon-nucleon collision scaling {i.e. Rqp ~ 0.33). This 
Pt value establishes the extent to which the centrality dependent enhancement 
of baryon production persists. 

The particle-type dependence of v 2 and Rep provides a stringent test for mod- 
els of heavy-ion collisions. In particular the larger values of A + A v 2 compared 
to their smaller suppression manifested in R C p suggests that for p T < 4.0 GeV/c 
a particle production mechanism beyond the framework of energy loss and frag- 
mentation exists in central Au+Au collisions. The particle- and p?- dependence of 
v 2 , and Rcp are consistent, however, with expectations based on the hadroniza- 
tion of a bulk partonic matter by coalescence or recombination. As such, the 
constituent- quark- number scaled v 2 reflects the anisotropy established in a par- 
tonic stage and provides strong evidence for the existence of a quark-gluon plasma 
in Au+Au collisions at RHIC. 
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CHAPTER 1 



Introduction to Relativistic Heavy-Ion 

Collisions 

By colliding heavy nuclei at relativistic energies scientists are able to test 
the nature of nuclear matter at high temperature and density, to produce condi- 
tions similar to those thought prevalent in the early universe, and to search for 
previously unstudied states of nuclear matter. In this chapter, we discuss the 
essential components of the theory thought to govern heavy-ion collisions, and 
we introduce the analysis topics that will be presented in this thesis. 

1.1 QCD— Asymptotic Freedom and Confinement 

Matter is made of leptons, quarks, and force mediators. Quarks, the building 
blocks of nucleons (and all hadronic matter), carry a property analogous to elec- 
tric charge called color. The theory that describes the forces between colored 
objects and that is thought to be the correct theory for strong interaction is 
called quantum chromodynamics (QCD). In QCD, just as the electromagnetic 
force is carried by photons, the color force (or strong force) is carried by gluons. 
However, whereas photons carry no electric charge, gluons do carry color charge 
so they can interact directly with each other, and whereas the electrodynamic 
coupling constant a = t|», the strong coupling constant a s can be larger than 
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one. As a consequence of the direct gluon-gluon coupling the effective coupling 
constant for the strong force becomes smaller at shorter distances. This effect 
is known as asymptotic freedom. Asymptotic freedom means the force between 
quarks is stronger at larger distances so quarks seem to remain confined to a 
small (~1 fm 3 ) region in colorless groups of two (mesons) or three (baryons). Be- 
cause the effective strong coupling is only small at short distances, perturbation 
theory can only be used with QCD for interactions involving large momentum 
transfers (i.e. hard processes). Although perturbative QCD (pQCD) is in very 
good agreement with experimental observations involving hard processes (see Fig- 
ure for example [UCOOj). it cannot be used to calculate QCD predictions for 
the processes that dominate the universe at present: soft processes 



Figure 1.1: Differential cross-sections for single jet production at pseudo-rapidity 
Tj = as a function of the jet transverse momentum pj- in proton (anti-)proton colli- 
sions. Jets are somewhat collimated sprays of particles produced when quarks or gluons 
collide, transfer (and carry away) a lot of momentum, and then fragment into a spray 
of hadrons. The curves represent pQCD calculations for the collisions at center-of-mass 
energy Vs = 630 and 1800 GeV. 
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Explicit QCD Lagrangian calculations of the force between quarks can only 
be made in the limits of weak and strong coupling. To understand the behavior 
of colored objects where pQCD is not a valid approximation, physicists rely 
on numerical path integrals of the QCD Lagrangian on a discretized lattice in 
four-dimensional Euclidean space-time. It is the formulation of Lattice QCD 
with a strong coupling approximation that first demonstrated how quarks are 
confined |Wil74j . 

In principle, the lattice formulation of QCD can be used to perform numerical 
calculations for all physical regimes. In practice, however, there are regimes 
where approximations used to simplify the calculations fail and the computations 
become technically very challenging. 

1.2 Deconfined Quark Matter 

In the strong coupling regime the energy required to separate two quarks increases 
linearly with the distance between them. As a result, we have never observed 
deconfined quarks: a deconfined quark is taken as one that can move in a volume 
much larger than the volume of a proton. Recent advances in the formulation of 
thermodynamical lattice QCD at finite temperature and density however, sug- 
gests that when sufficiently high temperature and density are reached, quarks 
become effectively deconfined. Figure 11.21 [Kar02 shows that the ratio of the 
energy density scaled by T 4 (where T is the system temperature) e/T 4 quickly 
increases at a critical temperature Tq- The magnitude of e/T 4 reflects the num- 
ber of degrees of freedom in the thermodynamic system. The rise corresponds to 
a transition in the system to a state where the quarks and gluons have become 
relevant degrees of freedom. 
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Figure 1.2: The energy density in QCD from lattice calculations. When the temper- 
ature T reaches the critical temperature T c , the number of degrees of freedom rapidly 
rises indicating that quarks and gluons become relevant degrees of freedom. The arrows 
represent the Stefan-Boltzmann values for asymptotically high temperature. 

The idea of a new state of matter where deconfined quarks and gluons are the 
relevant degrees of freedom is not new. In 1973, shortly after asymptotic freedom 
was shown to arise from QCD theory [GW73, Pol73j, deconfined quark matter 
was postulated as the true state of nuclear matter at high energy density at the 
center of neutron stars |CP75j : 

A neutron has radius of about 0.5-1 fm, and so has a density of about 
8 x 10 14 gm/cm 3 , whereas the central density of a neutron star can 
be as much as 10 16 — 10 17 gm/cm 3 . In this case, one must expect the 
hadrons to overlap, and their individuality to be confused. Therefore, 
we suggest that there is a phase change, and that nuclear matter at 
such high densities is a quark soup. 

Later, in the fall of 1974, at a workshop on heavy-ion collisions, T.D. Lee discussed 
the need for a physics program to study quark matter |Lee75j : 

Hitherto, in high-energy physics we have concentrated on experiments 
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in which we distribute a higher and higher amount of energy into a 
region with smaller and smaller dimensions. In order to study the 
question of "vacuum," we must turn to a different direction; we should 
investigate some "bulk" phenomena by distributing high energy over 
a relatively large volume. 
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Figure 1.3: The energy density e and pressure P, scaled by T 4 from a statistical 
model RL03j. The various lines show results for different hadron mass spectra. The 
results show an increase in the energy density degrees of freedom at a critical temper- 
ature near T = 158 MeV. 



Not all conceptualizations of the cross-over from hadronic degrees of freedom 
to a new form of matter relied on QCD or the knowledge of quarks. In 1951, 
Pomeranchuk postulated an upper limit to the temperature of hadronic mat- 
ter based on the finite size of hadrons |Pom51j . In the late sixties, Hagedorn's 
approach involving a self-similar hadronic resonance composition pointed to a 
similar limit |Hag65]. We now believe these limits reflect a transition to a state 
of matter with quarks and gluons as deconfined constituents. Figure 11.31 shows 
the scaled energy density e/T 4 and scaled pressure P/T 4 derived from a statistical 
model of a hadronic gas 
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1.3 Goals of Heavy-Ion Physics 



The creation and study of bulk matter made of deconfmed quarks and gluons (i.e. 
a quark-gluon plasma or QGP) was one of the prime motivations for building the 
Relativistic Heavy-Ion Collider (RHIC). The interaction of high-energy, colliding 
beams of heavy nuclei generates matter of extreme density and temperature. The 
temperatures and densities reached are expected to be similar to those thought 
to have prevailed in the very early universe, prior to the formation of protons 
and neutrons. The observation and study of matter in these conditions will 
be relevant to the nuclear physics community, the astrophysics community and 
the high-energy physics community. One also expects this research to have a 
significant impact on many in the general public since the nature of our universe 
at the earliest stages and the transitions that produced the matter we are familiar 
with today are interesting to most naturally curious or inquisitive people. 

By colliding large nuclei at high energy a window is opened onto an asymptotic 
regime of QCD. The exploration of this region of the QCD phase diagram is 
an exciting scientific endeavor. Many questions will be addressed in heavy-ion 
research programs: How well does the system thermalize? In the early universe, 
how did matter hadronize? Is there a first order phase transition, second order 
phase transition or smooth cross-over? How is fragmentation affected by the dense 
system created in the collisions? What is the role of chiral symmetry breaking in 
the transition from deconfmed partons to hadrons? The measurements presented 
here provide insight into how well the matter created at RHIC thermalizes and 
how it subsequently hadronizes. 

Learning about dense nuclear matter is also important to the astrophysics 
community. Heavy-ion physics can potentially provide insight into the structure 
of neutron stars (i.e. their mass-radius relationship, their thermal evolution, their 
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upper mass limit). In addition, reaching a better understanding of dense nuclear 
matter will help determine whether a new class of stars, quark stars, are likely 
or unlikely to exist in our universe. 

Perhaps the most exciting discoveries made will be those that are least ex- 
pected. The heavy-ion collisions at RHIC constitute an exploration into the 
unknown and one should be ready to be surprised. We do not know, for example, 
what, if any, exotic states may be produced in the hadronization of bulk quark 
matter. Candidates include multi-quark states, exotic atoms, and large droplets 
of strange-quark matter. 

1.4 Experimental Observations 



Figure 1.4: A sketch of the expected evolution of a relativistic heavy-ion collision. 

Figure ETT1 depicts the space-time evolution of a heavy-ion collision. Four pos- 
sible stages of the evolution are shown: a pre-equilibrium stage, an equilibrated- 
deconfined-parton stage, an interacting-hadron-gas stage and finally a free-hadrons 
stage. The experiments at RHIC detect hadrons in the free-hadron stage of the 
collision evolution. Probing the early stage of the collision evolution with parti- 
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cles measured in this final stage is a significant challenge. In this thesis we will 
present measurements thought to be sensitive to the early part of the collision 
evolution and to a possible deconfined-parton phase. 



1.4.1 Initial Conditions 

It is not known a-priori that an equilibrated-deconfined-parton phase can be cre- 
ated by colliding heavy ions in the laboratory. The large energy densities reached 
in central collisions (i.e. head-on collisions) however, significantly surpass esti- 
mates of the energy densities needed to reach the deconfinement phase transition. 
The initial energy density e of the produced medium can be determined using 
the Bjorken estimate |Bjo83| 



where (-^ L ) y=0 is the number of hadrons per unit rapidity produced at mid- 
rapidity Wh is the average energy of the hadrons, Ra is the nuclear radius, and r 
is the formation time of the medium. The formation time is not known but is gen- 
erally taken to be approximately one fm/c. The density of normal nuclear matter 
is approximately 0.16 GeV/fm 3 . Lattice calculations predict that the phase tran- 
sition to deconfined quarks and gluons occurs near 1.0 GeV/fm 3 . The Bjorken 
estimate for the initial energy density in central Pb+Pb collisions with y/s NN = 
17 GeV at the CERN-SPS experiment is 3.5 GeV/fm 3 Sat 03 . The estimate from 
RHIC for central Au+Au collisions at ^/s^ = 130 GeV is 4.6 GeV/fm 3 |Zaj02| . 
For the top RHIC energy (^s^ = 200 GeV), e ~ 5.0 GeV/fm 3 [DE03j. These 
estimates of e far exceed the energy density thought necessary to generate de- 
confined partonic matter. Given these large densities, collective behavior due to 




(1.1) 
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multiple interactions is expected. An important question to ask then is; are the 
interactions copious enough and rapid enough to thermalize the dynamic and 
expanding matter created in the laboratory? Answering this question will be a 
challenge to the experiments at RHIC. Figure 11.51 (left) shows how the rapidity 
density per participating nucleon pair increases as a function of y/s NN jBaz03| . 
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Figure 1.5: Left: Charged particle rapidity density at mid-rapidity scaled by the 
number of participating nucleon pairs 0.5 x N part versus y/s NN Raz03 . Right: Scaled 
pseudo-rapidity density for Au+Au and p+p collisions at y/s NN = 19.6, 130, and 
200 GeV versus N part |Bac03bj . 



In addition to the center-of-mass energy ■ s /s NN , the initial conditions of heavy- 
ion collisions also depend on the centrality of the collision. An off-axis nucleus- 
nucleus collision will have a smaller number of participating nucleons (N part ), a 
smaller system size, and a smaller initial energy density. Figure IT31 (right) shows 
the rapidity density per participating nucleon pair versus N part |Bac03b| . We 
also note that for nuclei colliding off-axis, the overlap region will be asymmetric. 
In Figure 11.61 we plot the overlap density for Au nuclei colliding with impact 
parameter b = 5 fm. A Woods-Saxon distribution is used for the density profile 
of the Au nuclei. 
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Figure 1.6: The overlap density for Au nuclei colliding off-axis. The beam directions 
are in and out of the of the page. The two large circles represent the outline of the 
incoming or outgoing Au nuclei. The impact parameter b — the distance between the 
center of the two colliding nuclei in the azimuthal plane — is 5 fm. The reaction plane 
is by the beam axis and the vector connecting the centers of the two nuclei. 

Most observables in heavy-ion collisions are integrated over the azimuthal 
angle and, as such, they are insensitive to the azimuthal asymmetry of the initial 
source. In this thesis we discuss measurements sensitive to the conversion of 
the initial spatial anisotropy to a final momentum-space anisotropy. The spatial 
anisotropy can be quantified by estimating the eccentricity e of the initial source, 



Extracting the mean eccentricity of the initial source (e) for a given centrality 
interval is helpful for understanding the event-by-event anisotropy in the final 
state momentum distributions. Analytic of initial eccentricities can be found in 
Appendix 1X1 



e = 



{y 2 - x 2 ) 

(y 2 + x 2 ) 
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1.4.2 Event-by-event Momentum-space Anisotropy 

Anisotropy in the distribution of a particle in momentum-space is thought to be 
sensitive to the early stage of the collision system. The anisotropy of the source 
will be largest immediately after the collision occurs. As the system evolves, 
the spatial anisotropy is converted by multiple interactions into a momentum- 
space anisotropy. With time, the interactions will cause the spatial distribution 
to become more isotropic. For this reason, it's believed that the final azimuthal 
momentum-space anisotropy is primarily built up in the initial moments of the 
system's evolution. Figure 11.71 shows the evolution of the source shape calcu- 
lated from a model where the collision system is described by hydrodynamic 
equations [KSHOO . 

The azimuthal anisotropy of the transverse momentum distribution for a par- 
ticle can be described by expanding the azimuthal component of the particle's 
momentum distribution in a Fourier series, 



The harmonic coefficients, v a , are anisotropy parameters, px, y, and are the 
respective transverse momentum, rapidity, and azimuthal angle for the particle, 
and tyRp is the reaction plane angle [PV98] 1 . The second coefficient v 2 (custom- 
arily called elliptic flow) measures the elliptic component of the anisotropy. Due 
to the shape of the source created in an off-axis collision, v-i is the largest and 
most studied of the anisotropy parameters. Figure 11.81 shows the energy depen- 



lr The reaction plane is defined by the beam axis and the vector connecting the centers of 
the two colliding nuclei. For high energy collisions, in the laboratory reference frame the Au 
nuclei are Lorentz-contracted along the beam axis. As such, the vector connecting the colliding 
nuclei is nearly perpendicular to the beam axis and the reaction plane can be characterized by 
its azimuthal angle. 
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Figure 1.7: The evolution of the source shape is shown from a model where a heavy-ion 
collision is treated as a hydrodynamic system [KSHflOj . The initial shape is extended 
out-of-plane. By 8 fm/c after the formation time (r — To), the shape has deformed to an 
in-plane extended source. In this model, the anisotropy in momentum-space measured 
by V2 is dominated by the early stages. 
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Figure 1.8: The integrated anisotropy parameter V2 near mid-rapidity for mid-central 
events (roughly 12-34% central) plotted versus collision energy |Alt03j . 

dence of t>2 for charged particles near mid-rapidity. For this energy range t>2 is 
positive and rising monotonically with the nucleon-nucleon center of mass energy. 
At lower energies v 2 is neg ative H3a02 . 

Multiple interactions are necessary to develop a momentum-space anisotropy 
from a coordinate-space anisotropy. If each nucleon-nucleon collision is indepen- 
dent, the final momentum distribution will represent a superposition of random 
collisions and will therefore be isotropic. The azimuthal momentum-space distri- 
bution of charged hadrons with 2.0 < px < 6.0 GeV/c, for three centrality inter- 
vals, is shown in Figure EH (left). Figure EH (right) shows how V2 for charged 
hadrons changes with px (differential ^2). The magnitude of t>2 is smallest in 
central events because the initial eccentricity e is smaller. 

The large saturated values of v 2 at high p T are a surprising result from RHIC. 
Although hydrodynamic models predict a monotonic increase of differential v%, it 
is believed that hydrodynamic models must fail at higher values of pt where their 
assumptions become invalid. The measurement of a large t>2 at high p? gives rise 
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to the question: "how does the initial spatial anisotropy manifest itself in the 
distribution of high pt particles?" One explanation is that high energy partons 
lose energy as they pass through the matter created in the collisions. Since the 
source is asymmetric, the amount of energy loss will depend on the direction the 
parton travels. As such, energy loss can lead to a momentum- space anisotropy 
that reflects the initial spatial anisotropy of the source. We will discuss energy 
loss and the suppression of high p T particle production in section 11.4.31 
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Figure 1.9: Left: The distribution of charged particles with 2.0 < pr < 6.0 GeV/c in 
the azimuthal plane from Au+Au collisions at y/s NN = 130 GeV. The 0-10%, 10-31%, 
and 31-77% represent different classes of centrality where 0-10% is the most central. 
Right: The differential V2 for in three centrality intervals [Adl03a . 



Figure IT. 1 01 shows the differential v?, at mid rapidity for identified particles at 
low pt {pt < 1 GeV/c) where particles can be identified by their energy loss in the 
detector gas jAdl OlJ. The hydrodynamic models predict a mass-ordering for ellip- 
tic flow with less massive particles having larger elliptic flow for all values of pt- 
The large and its mass-ordering at low pt are consistent with the hydrodynamic 
limit for the conversion of spatial anisotropy to momentum anisotropy (where lo- 
cal thermal equilibrium has been assumed) jHKHOll IOU921 ISor99l ITLSOlj . At 
intermediate pt (1.5 < pr < 4 GeV/c), however, while hydrodynamic models 
predict a monotonic increase, the charged hadron t> 2 saturates at a value ap- 
proximately independent of px- After the first year of RHIC data taking, the 
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Figure 1.10: The differential elliptic flow for identified particles at mid-rapidity from 
Au+Au collisions at ^ys NN = 130 GeV {AdlOlj . The curves represent fits to hydrody- 
namic inspired parameterizations. 

particle-type dependence of v-i in the high momentum region remained an open 
question. In this thesis we present measurements of for Kg and A + A at 
mid-rapidity from Au+Au collisions at y/s NN = 200 GeV/c that extend up to 
p T ~ 6 GeV/c. 

1.4.3 Nuclear Modification of Particle Production 

Like t>2, high p? hadron production — presumably through scatterings of partons 
involving large momentum transfer — is also thought to probe the early stage 
of heavy-ion collisions. High-energy partons passing through dense matter are 
predicted to lose energy by induced gluon radiation jGP90l IBSZ001 IGVW 03J. 
Since the total energy loss depends on the color charge density of the medium, 
nuclear modification of the high p? particle yields can probe the dense, perhaps 
deconfmed-partonic matter created by the collision. 

Partonic energy loss or jet-quenching can be studied by measuring the modi- 
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fication of particle production in nuclear collisions. A nuclear modification factor 
can be formed by taking the ratio of the particle yields in nucleus-nucleus colli- 
sions and the particle yields in proton-proton collisions. The ratio is then scaled 
by Taa = (Nbmar y )/o"i^ to account for the trivial increase in the yield with the 
system size: 



Raa{P' 



T 



d 2 n AA / dp T dr] 



;i.4) 



T AA d 2 a NN /dp T dr]' 

where rj is the pseudo-rapidity and Nonary is the number of binary nucleon-nucleon 
collisions. In the absence of nuclear effects, at high p?, Raa is expected to be 
unity. In the low pt region the yield is not expected to scale with N binary . The 
p<r-scale where the high p? regime begins is an experimental observable that our 
measurements will address. 
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Figure 1.11: The nuclear modification parameter Raa for charged hadrons |Ada 03b 
and for neutral pions |Adl03f] in Au+Au collisions at y/s NN = 200 GeV. The panels 
show different Au+Au collision centrality intervals with the most central in the top left. 
The curves show model calculations based on partonic energy loss and other nuclear 
effects [WftTi()3j . 



Figure IT. 1 II shows Raa for charged hadrons and neutral pions 7r° from Au+Au 
collisions at y/s NN = 200 GeV. The high yields in central Au+Au collisions are 
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suppressed with respect to Nonary scaling. The suppression is largest for central 
collisions while the yields in peripheral collisions are consistent with expectations 
from Nbinary scaling. The suppression is approximately independent of pt for 
Pt > 3 GeV/c for tt° and for px > 5 GeV/c for charged hadrons. 
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Figure 1.12: Nuclear modification of charged particle production measured at 
y/s NN = 200 GeV from central Au+Au events (0-5% of the collision cross-section) 
compared to peripheral Au+Au events (40-60% and 60-80% of the collision 



cross-section) |Ada 03b . 



Like Raa, the ratio of the yields in central and the yields in peripheral colli- 
sions (Rcp) also can measure nuclear modifications to particle production: 



central 



R C p{p 



T 



[dn/ (N binary <ipr)] 
[dn/ (N binary <ipT)] peripheral 



;i.5) 



R central_ 



When Raa for peripheral events follows N binary scaling, R C p 
ratio Rcp typically has smaller systematic uncertainties than Raa and does not 
require the measurement of a p+p reference spectrum. The charged hadron Rcp 
in Figure IT. 1 21 shows a suppression of particle yields in central events compared to 
scaled peripheral events. The suppression of charged hadrons is roughly constant 
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for pt > 5 GeV/c. The dependence on particle-type of the suppression and the 
Pr-scale for its onset remained an open question after the first year of RHIC 
collisions. In this thesis we present the measurement of Rqp for Kg and A + A 
from Au+Au collisions at y/s NN = 200 GeV up to pr ~ 6 GeV/c. 

As mentioned in section H.4.2[ energy loss can also manifest itself in Vi- By 
suppressing the yield of large px particles more in the out-of-plane direction than 
the in-plane direction 2 , energy loss can cause an anisotropy in the final momentum 
distribution. The particle-type dependence of t>2 and Rqp will be a powerful test 
of the energy loss hypothesis. If energy loss governs the development of v 2 and 
Rcp then we expect either no particle- type dependence or we expect the particle 
with the larger v-i to also have a larger suppression. 



1.4.4 Other Observations 




Number of participants Numberof participants 



Figure 1.13: Left: Centrality dependence of the percolation scale Qs where char- 
monium suppression is thought to set in for RHIC energy and SPS energy |DF P02 . 
Right: The ratio of the measured and expected J/tj) yield showing a step like suppres- 
sion pattern as the number of participants in the collisions system increases |Sat03j . 



Other important observations made in heavy-ion collisions include the sup- 



2 The in-plane and out-of-plane directions are perpendicular to the beam axis. The in-plane 
direction lies along the vector connecting the colliding nuclei. 
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pression ofJ/ip production, strangeness enhancement, the coincidence of particle 
ratios with statistical model predictions, the enhancement of baryon production, 
and the reduction of the net baryon number. Figure IT. 131 (right) shows the ratio 
of the expected and measured J/^ yields versus N part . The step like behavior was 
interpreted as being caused by the dissolution of successive charmonium states in 
a new form of matter created in heavy-ion collisions. First the ip' and xc dissolve 
and then the J/ip dissolves. Figure IT. 131 fleft) shows the percolation scale where 
the dissolution of different charmonium states should set in versus N part [DFP02J. 
Figure 11.141 shows the enhancement of the (anti-)proton to pion ratio at inter- 
mediate pt in central Au+Au collisions relative to e + + e~, p+p, or peripheral 
Au+Au collisions. The enhancement of baryon production will be studied further 
in this thesis. Figure 11.151 shows the enhancement of strange particle production 
in heavy-ion collisions relative to p+Be collisions. The enhancement increases 
with the strange quark content; i.e. f2~(sss) > H~(dss) > A(uds) |Fan02| . 
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Figure 1.14: The proton/pion and anti-proton/pion ratio |Adl03ej . 
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Figure 1.15: Strange particle yields scaled by N part normalized by p-Be colli- 
sions |Fan02| . 



1.5 Thesis Outline 



In this thesis, we present the first measurement of Rcp for K s and A + A for 
Au+Au collisions at y/s NN = 200 GeV and the first measurements of v 2 for K s 
and A + A for Au+Au collisions at y/s NN = 130 and 200 GeV. Our emphasis 
is on probing the early stage of heavy-ion collisions, mapping out the transition 
between pt regions (i.e. soft, intermediate, hard, etc.), and understanding how 
hadronization modifies the observables we measure. In mapping out the pt re- 
gions we hope to learn what processes dominate particle production within each 
region. Studying the variation in yields with centrality (Rcp) and azimuthal 
angle (v 2) for different particle species will help us understand the hadronization 
mechanisms in heavy-ion collisions. This information will be helpful for charac- 
terizing the matter created in heavy-ion collisions. 

In Chapter 2 we will discuss the facilities used to study heavy-ion collisions. 
The Relativistic Heavy-Ion Collider (RHIC) will be described, an introduction 
to particle tracking detectors will be given, and the Solenoidal Tracker at RHIC 
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(STAR) detector system will be reviewed. Chapter 3 contains details of the 
analysis methods. In Chapter 4 we present the results of the analysis and in 
Chapter 5 we discuss these results, draw conclusions, and present an outlook for 
future work. In the appendices we include a description of the coordinates system 
in the transverse plane, calculations of the nuclear overlap density for Au+Au 
collisions, definitions for the kinematic variables used in this thesis and a list of 
STAR collaborators and institutions. 
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CHAPTER 2 



Experimental Set-up 
2.1 The Relativistic Heavy-Ion Collider 

The Relativistic Heavy-Ion collider (RHIC) at Brookhaven National Lab 
(BNL) is designed to collide counter-rotating heavy-ion beans at energies up 
to 100 GeV/u. RHIC is the first facility to collide heavy-ion beams. The center- 
of-mass energy for these collisions is roughly a factor of ten times greater than the 
highest energies reached with the previous fixed target heavy-ion experiments. 
Parameters for existing and future relativistic heavy-ion facilities are given in 
Table 12.11 RHIC consists of two concentric rings of super-conducting magnets 
(cooled to below 4.6 degrees Kelvin) that focus and guide the beams and a ra- 
dio frequency (rf) system that captures, accelerates and stores the beams. The 
ring's diameters are approximately 1.22 km. 

Figure 12.11 shows the BNL accelerator complex including the accelerators used 
to bring the gold ions up to RHIC injection energy. In the first of the Tandem 
Van de Graaff accelerators, gold ions in a charge state Q = — le accelerate to 15 
MeV. The ions then pass through a stripping foil (located between the Van de 
Graaffs) where electrons are knocked off so that their most probable charge state 
becomes Q = +12e. With their charge changed from negative to positive, the 
ions gain another 1 MeV/u of energy as they accelerate through the second Van 
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Figure 2.1: A diagram of the Brookhaven National Laboratory collider complex in- 
cluding the accelerators that bring the nuclear ions up to RHIC injection energy (10.8 
GeV/u for 197 Au 79 ). 
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Table 2.1: RHIC compared to existing and future facilities; A max is the maximum 
species mass number, Ep aa; is the maximum (equivalent) fixed-target beam energy 
per nucleon, y/s NN is the maximum center of mass energy per nucleon, y/s AA is 
the total center of mass energy, and Ay/2 is the rapidity gap from the beam to 
mid-rapidity |LR02| . 

de Graff, back to ground potential. On exiting the Tandem, the ions pass through 
a second stripping foil bringing their most probable charge to Q = +32e. They 
are then injected into the Booster synchrotron and accelerated to 95 MeV/u. A 
stripper foil in the transfer line between the booster and the Alternating Gradient 
synchrotron (AGS) increases their charge state to Q = +77e. In the AGS the 
ions are accelerated to 10.8 GeV/u. They are extracted from the AGS and passed 
through one final stripper foil where the remaining K-shell electrons are removed 
(Q = +79e). Finally, they are injected into RHIC where they are accelerated 
to top energy and can be stored for up to 10 hours. Table 12.21 lists important 
parameters for RHIC. 

Top Au+Au v /s^7 200 GeV 
Ave. luminosity C (10 hour store) ~ 2 x 10 26 cm~ 2 s _1 

Bunches per ring 60 

Gold ions per bunch 10 9 

Crossing points 6 

Beam lifetime (store length) ~ 10 hours 

RHIC circumference 3833.845 m 

Table 2.2: Nominal RHIC parameters for Au+Au collisions. 
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2.2 RHIC Experimental Program 



To date, RHIC has generated collisions between gold nuclei at y/s NN = 22, 56, 
130, and 200 GeV, between protons at y/s NN = 200 GeV, and between gold and 
deuterium nuclei at y/s NN = 200 GeV. Table 12.31 shows the luminosity achieved 
at the end of RHIC Run-2 (Au+Au collisions at ^/s^ = 200 GeV). The STAR 
experiment recorded integrated luminosities £ ~ 2.8 b" 1 and £ ~ 80 fj,b~ l for 
RHIC Run-1 (Au+Au collisions at y/s NN = 130 GeV) and RHIC Run-2 respec- 
tively. Most of the integrated luminosity comes late in the runs, after the collider 
is tuned. 

Bunches Ions/Bunch £ pea k £ave (store) Integrated 

[cm-V 1 ] [cm-V 1 ] £ [(fib)- 1 ] 

55 6 x 10 8 3.7 x 10 26 3.7 x 10 26 ~ 80 

Table 2.3: The performance of RHIC during the 2001 Au+Au run (Run-2). 

There are four experimental collaborations at RHIC; the PHOBOS collabo- 
ration with 107 members from 8 institutions, the BRAHMS collaboration with 
51 members from 14 institutions, the PHENIX collaboration with 328 members 
from 52 institutions, and the STAR collaboration with 293 members from 39 
institutions 1 . 

In this thesis we present an analysis of Au+Au collisions recorded by the 
STAR detector during the summer of 2000 and the winter of 2001. Approximately 
5 x 10 5 and 5 x 10 6 usable events were recorded at y/s NN = 130 GeV and ^s NN = 
200 GeV respectively. 



1 These numbers are taken from each collaborations author list as of July 2003 and do not 
represent the total number of people working on the experiment. There are, for example, over 
450 scientist and engineers working on the PHENIX experiment. 
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PHENIX event display 




STAR event display 



Two muon spectrometers cover the 
pseudo-rapidity region 1.1 < < 
2.4 and azimuth angle < <fi < 2n. A 
central spectrometer with two arms 
and tracking sub-systems (each sub- 
tending n / 2 radians) covers \rj\ < 
0.35. With a smaller acceptance and 
faster detectors the emphasis is on 
triggering on rarer probes, hadron 
identification and electron identifica- 
tion. 



A large acceptance solenoidal track- 
ing detector with particle identifica- 
tion covers the full azimuth (\<f>\ < tt), 
\rj\ < 2.0 and 2.5 < \rj\ < 4.0. 
Subsystems include a central TPC, 
two forward TPCs, a silicon vertex 
tracker and a barrel electromagnetic 
calorimeter. The emphasis is on 
global event characterization, reso- 
nance identification, fluctuations and 
event-by-event variables. 





PHOBOS event display 



BRAHMS detector 



Measurements of charged particles 
are made across a full solid-angle 
with a multiplicity detector. Two 
small acceptance spectrometer arms 
allow for particle identification at 
mid-rapidity. Multiplicity measure- 
ments across a broad range of rj and 
p T are emphasized. 



Designed to provide good particle 
identification across a broad rapidity 
and pt range (0 < y < 4; 0.2 < 
Pt < 3.0 GeV/c) with two small 
solid-angle spectrometers. Measur- 
ing particle production at forward 
angles is emphasized. 



Table 2.4: Summary of RHIC experiments. 
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2.3 Particle Tracking Detectors 



The primary detector used for the analysis presented in this thesis is the STAR 
Time Projection Chamber (TPC). The TPC is designed to do particle tracking 
which facilitates the identification of secondary vertices from weak decays (e.g. 
Kg — * 7r + + In the following we give an introduction to high-energy particle 
tracking technology [SKN03 . 

2.3.1 History of Particle Tracking 

Since the beginning of particle physics, when J.J. Thomson realized that the 
cathode rays he was studying were not "rays" but streams of subatomic charged 
particles instead, our understanding of the subatomic realm and the mechanics 
that governs it has depended strongly on our ability to detect the tracks of charged 
particles. Thomson was able to surmise the existence of electrons, measure their 
charge to mass ratio and even measure their velocity as they were emitted by a 
hot filament because he could see their trajectory as they passed through crossed 
electric and magnetic fields. During the years since Thomson's experiments in 
1897 many techniques have been developed to detect or visualize charged particle 
tracks — nuclear emulsions, cloud chambers, bubble chambers, spark chambers, 
streamer chambers, other gas detectors, solid-state detectors and so on. All of 
these techniques rely on very fast charged particles ionizing atoms as they pass 
through matter. The ionization left along the paths of the high-energy particles 
can then act as catalysts for reactions that leave an observable trace — such as, a 
bubble, a spark, condensation, or a charge cascade. 

Many of these particle detecting techniques have — because of their inherent 
limitations — been abandoned for the most part in favor of gas or solid state de- 
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tectors. Experiments at the newest colliders like RHIC or the LHC rely almost 
exclusively on these two techniques because they lend themselves well to trig- 
gering, high event rates and the digitization of huge amounts of data. Bubble 
chambers however, are of particular historic importance and produced a wealth of 
information from their inception in 1952 to well into the 1970s. Their importance 
to high-energy physics was acknowledged with a Nobel prize in 1960 awarded to 
Donald Glaser. Glaser struck upon the idea of a bubble chamber when he saw 
the tracks created by bubbles in beer. 

2.3.2 Bubble Chambers — Three Decades of Physics and Two Nobel 
Prizes 

Bubble chambers initially used liquid in a super-heated state to detect the 
ionization left along the tracks of high-energy charged particles passing through 
the liquid. In 1952 Glaser used diethyl-ether heated to ~ 100°C above its boiling 
point to build the first bubble chamber. The super-heated liquid, when struck by 
cosmic rays, began boiling violently and a photograph made using a fast camera 
showed tracks left behind by the high-energy charged particles created by a cosmic 
ray. It is presumed that after a high-energy particle causes the initial ionization 
along its path, heat generated by recombination is responsible for the boiling and 
bubble formation in the liquid. 

Improvements to this technique, including the use of pistons to create a sud- 
den pressure drop to induce bubble formation, led to greater precision and larger 
chambers that could be placed in a magnetic field. Perhaps the most notable im- 
provement to the bubble chamber came when Luis Alvarez substituted hydrogen 
for the ether used by Glaser. Alvarez's chamber produced much clearer tracks 
and this technological advance was of such value that he won the Nobel prize for 
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physics in 1968 for his work, the second Nobel prize awarded for work related to 
the development of the bubble chamber. 

In bubble chambers, the fluids in the chamber act both as the target and as 
the detector so different fluids — some cryogenic and some room temperature — 
were eventually used to suit the purpose of the experiment. Cryogenic liquids 
consisted of the simplest nuclei like, H 2 , D 2 , He, Ne, Ar and Xe while the room 
temperature "heavy liquids" like propane (C 3 H 8 ) and Freon (CF 2 C1 2 or CF 3 Br) 
offered short interaction lengths. The typical size of a bubble in a bubble chamber 
is ~ 10 (Km and the bubble density can be used to determine (3 = v/c for the 
passing particle. 

The advantages bubble chambers offered kept them in widespread use for 
three decades, from the 1950's to well into the 1970's. They had good spatial 
resolution (10 - 150 /iin), a large sensitive volume, 4n geometrical acceptance; and 
they permitted the use of a variety of materials as targets. Eventually however, as 
physics began requiring more complex triggers and as large- volume high-precision 
detectors demanding electronic data recording came in use, the bubble chambers 
disadvantages rendered them obsolete. The analysis of photographs was a tedious 
task requiring expensive projectors for scanning the images and the whole process 
was only modestly scalable so that only limited statistics could be achieved. 
Bubble chambers were also complicated to operate, required cryogenics, and were 
a safety concern. In addition, bubble chambers weren't compatible with particle 
colliders — the now dominant high-energy accelerator, they provided no triggering 
for low cross-sections and they had a relatively long sensitive time (~ 1 ms) which 
necessitates a lower beam luminosity. 
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2.3.3 Streamer chambers — a Precursor to Modern Gas Detectors 

The streamer chamber developed by G.E. Chikovani in 1963 — an improvement 
on the spark chamber — overcame some of the limitations of the bubble chamber 
and was the predecessor of the gaseous detectors of today. Like the bubble cham- 
ber however, the streamer chamber also relied on photographic film to record the 
tracks of streamers, placing a limit on the statistics available for analysis. 

The spark chamber uses a large potential across two parallel planes of elec- 
trodes to induce electrical breakdown — a spark — in gas between the electrodes. 
The ionization left from the passing of a high-energy particle acts as the cata- 
lyst for the spark. The spark chamber however, can only measure the position 
of the track in the direction parallel to the electric field to within the spacing 
of the electrodes. The streamer chamber overcomes this limitation by apply- 
ing a high- voltage pulse for a short duration (~ 15 ns). The strong electric field 
(~ 20 kV/cm) from the high- voltage pulse induces an incomplete spark discharge. 
These electron avalanches or streamers form all along the particles path and the 
radiation of the gas in the streamer plasma can be recorded optically. Streamer 
chambers were built with sensitive volumes of several cubic meters that recorded 
particle tracks in any direction with equal efficiency The density of the streamers 
can be used for particle identification up to particle momenta of ~ 1 GeV/c. 

The two major advantages of the streamer chamber over the bubble chamber 
are its ability to be triggered by external devices and its very short sensitive time 
(~ 1 fis). Eventually however, its use of photographic film, its limited spatial 
resolution (> 300 /jm) and its relatively long dead time (~ 300 ms) turned out to 
favor the gaseous detectors that would rely on electronic, not optical, recording 
techniques. 
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2.3.4 Today's Tracking Detectors 

Almost all tracking detectors, other than the solid state and gaseous detectors 
using electronic recording techniques have been abandoned. The modern-era 
detectors have shorter sensitive times and shorter dead times so that the beam 
intensity of the particle accelerator can be increased and greater statistics can 
be recorded. These newer detectors also tend to be easier to operate and have 
greater spatial resolution. 

2.3.5 Gas Detectors 

Most gas detectors — multi-wire proportional chambers (MWPC), drift chambers, 
straw tubes, cathode strip or pad chambers, time projection chambers (TPC) 
and micro-strip gas chambers (MSGC) — use the proportional counting mode of 
operation. In this mode, the electrons from the primary electron-ion pairs created 
by the high-energy charged particle, are directed in an electrostatic field toward 
a very high field region (10-100 kV/cm) surrounding an anode wire of small 
radius. In this region, the fast electrons gain enough energy to create secondary 
electron-ion pairs. Each new electron produced by ionization, in turn, creates 
more electrons-ion pairs; the development of this avalanche or cascade is called 
gas multiplication. Most of the electrons in the avalanche are created very close 
to the wire so they are collected within a few nanoseconds. The heavier ions— 
also predominantly produced near the wire — move more slowly across a larger 
potential difference. As they do so, they induce a signal that can be detected 
with an amplifier and used for position and energy loss (dE/dx) measurements. 

The mode of operation of a gaseous detector is determined by the response of 
the ions and gas to the field strength surrounding the anode. In the proportional 
mode of operation the field strength is great enough to induce gas amplification— 
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typically 10 4 -10 6 times — but is not so great that it leads to complete breakdown 
or non- negligible space charge effects: caused by the build up of longer-lived posi- 
tive ions. In the proportional mode of operation the signal is proportional to the 
number of primary electron-ion pairs which is in turn proportional to the energy 
lost by the traversing particle. The measured dE/dx can then be used for particle 
identification (PID). 

A multi-wire proportional chamber consists of planes of independent 
wires — typically spaced 1-2 mm apart — set between two planes of cathodes at 
a distance of 3-4 times the wire spacing. A negative voltage is applied to the 
cathodes and the wires are held at ground. Each wire then acts as a proportional 
counter for primary electrons-ion pairs left along a particles track. The distance 
from cathode to anode is typically about 1 cm and the wire diameter should be 
about 20-50 /zm. The spatial resolution is given by dj y/l2 = 300 — 600 jum, 
where d is the wire pitch or spacing. For this invention G. Charpak was awarded 
the 1992 Nobel prize in physics. 

A drift chamber is a multi-wire proportional chamber with a large wire 
pitch — from several centimeters up to 50 cm but more typically 5 cm. Track 
position is determined by measuring the time electrons need to reach the anode 
wires. The speed of the electron drift depends on the gas used and the pressure 
in the chamber and is typically ~ 5 cm/jus so that a timing resolution of 1 ns 
gives a spatial resolution of ~ 50 /im. Different geometries and configurations 
can be used in order to create constant fields pointing toward the anode wires. A 
straw tube is a drift chamber composed of an individual straw shaped cathode 
(diameter of ~ 5 mm) with a single anode wire in the center. A "continuous" 
tracker can be constructed by packing many layers of straw tubes together. Straw 
tube detectors tolerate high loads because they don't use a common gas volume 
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and can achieve a resolution of about 150 jum with coarse time measurements. 

A time projection chamber (TPC) is a drastic variation on a simple drift 
chamber. A TPC consist of a large three-dimensional gas filled vessel with readout 
detectors on a wall at the end of a drift volume. The readout detectors are 
usually cathode pad chambers. A strong electric field across the TPC produced 
by a cathode on the wall opposite the readout planes creates the drift field. When 
a charged particle creates electron-ion pairs within the drift volume the strong 
drift field prevents them from recombining. The much lighter electrons move 
quickly toward the readout chambers. The drift field is chosen so that it is not 
strong enough to create secondary electron-ion pairs: typically hundreds of volts 
per centimeter. The readout chamber is separated from the drift volume by a 
gating grid. The gating grid is a plane of wires that electrostatically separates 
the amplification region from the drift region. The gating grid prevents the ions 
created in the amplification region from getting back into the drift region and 
allows for triggering of the detector; when an interesting event occurs the gating 
grid wires are set to voltages that allow electrons to pass through. 

The TPC readout chambers typically consist of an anode wire plane between 
a ground wire plane and a cathode pad plane. The signal induced on the anode 
wires is typically detected via image charges on several nearby pads. The posi- 
tion of the electron-ion cascade in the anode wire direction can be determined 
precisely by fitting a modified Gaussian to the signals on several consecutive 
pads. This measurement gives two transverse coordinates and the drift time 
gives the third coordinate, making the TPC a fully three-dimensional detector. 
Unlike other historic three-dimensional detectors however, such as the bubble or 
streamer chambers, the TPC is read completely electronically. The TPC also has 
the advantage that it has no pulsed very high-voltages and is fast compared to 
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historic detectors: its speed is determined by the maximum drift time which, for 
large chambers is ~ 100 /is. 

Inhomogeneities in the drift field and effects due to magnetic fields however, 
can distort the drift path of the electrons and further degrade the resolution. 
The electron clouds also diffuse at a rate of hundreds of /xm/ y / cm due to elastic 
rescattering in the gas as they drift toward the readout chamber. The TPC 
requires careful tuning of the drift field and a high degree of gas purity. Many 
parameters like drift length, track angle, or the number of primary ions affect the 
spatial resolution but a typical value is ~ 500 /xm. 

2.3.6 Solid State Detectors 

Solid state detectors — silicon micro-strip detectors, silicon pixel detectors and 
silicon drift detectors — offer very good resolutions of ± 10 - 100 /xm and are now 
in common use. Every detector planned at the Large Hadron Collider at CERN 
(LHC) will use trackers based on silicon devices. Silicon detectors require only 
3.6 eV of energy from the traversing particle to create an electron-hole pair. That 
is roughly one order-of- magnitude less than gas detectors require (~ 30 eV). This, 
along with silicon's higher density, means that the number of electron-hole pairs 
created by a minimum ionizing particle (MIP) in silicon is much greater than the 
number of electron-hole pairs created over the same distance in gas. In 1 jum of 
silicon a MIP produces ~ 100 charge pairs. To produce that much charge in gas 
would require several centimeters. As a result, unlike gaseous detectors, silicon 
detectors don't require signal amplification inside the detector; in a typical silicon 
detector a MIP will produce 20 - 30 thousand electrons. 

Typically, silicon detectors are built using ~ 300 /tm thick, high resistivity, 
n-doped silicon plates with a thin p-doped layer on one side. A reverse bias 
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voltage — positive on the n-side and negative on the p-side — is applied to deplete 
the silicon of free charge carriers and to create an electric field that will cause 
the electrons and holes to drift to opposite surfaces where readout structures 
are organized. The highly developed state of silicon technology allows for the 
production of many different readout structures. 

The readout structures for both silicon strip and silicon pixel detectors are 
layers of aluminum applied to the surface of the silicon. Silicon strip detectors 
use a solid layer of aluminum on the n-doped side of the silicon and a sequence 
of aluminum strips on the p-doped side. The strips typically have a pitch (d) 
of ~ 50 fim. The resolution for this pitch is ~ 15 fim or d/y/l2. The charge 
collected on the strip is electronically integrated and read out as an analogue or 
digital signal. If strips are placed on both sides of the silicon — a double sided 
silicon strip detector — two coordinates can be measured simultaneously. 

The silicon pixel detector uses pixels instead of strips and has the advantage 
that it is a true two dimensional micro-detector. Amplifier circuitry however, 
needs to be connected to each pixel which typically has a surface area of only 
~ 50 x 400/im 2 . This is done using specially designed readout chips that are 
bump-bonded to the detector silicon. Silicon pixel detectors are fast, have very 
low noise (small capacitance) and have excellent pattern recognition for high 
particle densities. They require however, a large number of readout channels 
(100 million is not uncommon), are very fragile and offer many technological 
challenges in development. 

Silicon drift detectors have two-dimensional capabilities but by their design 
avoid the large number of channels required by silicon pixel detectors. Silicon drift 
detectors use a silicon wafer, with an array of anodes arranged at one edge and 
cathodes at the other. An electric field drifts the primary electrons — from the 
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track of a passing particle — through the silicon, toward the array of anodes. A 
typical drift speed is ~ 15 mm//xs. The anode position along the edge of the 
wafer and the drift time give two coordinates for the position of the track. The 
third coordinate is given by the position of the wafer and, like a spark chamber, 
is only known to within the thickness of the wafer. The signal on the anodes can 
be read out at ~ 40 MHz — a very high frequency — but the time for the electrons 
to drift to the anodes (~ 5 /is), makes it a relatively slow detector. In addition, 
these detectors require very precise climate control because of the dependence 
of the drift time on temperature. The high resolution of all these solid state 
detectors however, makes them ideal for constructing vertex chambers that are 
particularly useful for detecting heavy-flavor particles. 

2.4 The STAR Detector System 




Figure 2.2: View of the STAR detector system. 
The STAR detector jAck03) (Figure 12. 2j) is an azimuthally symmetric, large 
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acceptance, solenoidal detector designed to measure many observables simulta- 
neously. The detector consists of several subsystems and a large Time Projection 
Chamber (TPC) located in a 0.5 Tesla solenoidal analyzing magnet. 
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Figure 2.3: Cutaway of the STAR detector in its 2001 configuration; including a par- 
tial installation of the electromagnetic calorimeter (EMC), the temporary ring-imaging 
Cherenkov detector (RICH), and a time-of-flight detector (ToF) prototype. 



The layout of the STAR detector system as it was for Run-2 is shown in 
Figure 12.31 The active subsystems included two RHIC-standard zero-degree 
calorimeters (ZDCs) that detect spectator neutrons, a central trigger barrel 
(CTB) that measures event multiplicity, a ring-imaging Cherenkov and time- 
of-flight detector that extend particle identification to higher p?, 10 percent of 
the full barrel electromagnetic calorimeter to measure photons, electrons and the 
transverse energy of events, and four tracking detectors. The tracking detectors 
are the main TPC, two forward TPCs, and the silicon vertex tracker (SVT). 

The TPC is STAR's primary detector |And03j and can track up to ~ 4 x 10 3 
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particles per event. For collisions in its center, the TPC covers the pseudo- 
rapidity region \rj\ < 1.8. It can measure particle px within the approximate 
range 0.07 < px < 30 GeV/c. The momentum resolution 5p/p depends on rj 
and pt but for most tracks Sp/p ~ 0.02. The full azimuthal coverage of the 
STAR detector (— n < (f> < n) makes it ideal for detecting weak decay vertices, 
reconstructing resonances and measuring w 2 and other variables requiring event- 
by-event characterization. 
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Figure 2.4: The pt reach of STAR's particle identification capabilities with the 2001 
detector setup. 



Figure l2~4l illustrates the STAR detector's particle identification capabilities 
during Run-2. These capabilities will be further enhanced with detector upgrades, 
larger data samples, and more advanced triggering to select rare events. Most 
of the measurements illustrated in Figure 12.41 are limited in px coverage by the 
statistics available. Using the topology of their weak decays in the TPC, the Kg 
and A(A) were identified across the largest pr range (0.3 < px < 7.0 GeV/c). The 
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kinematic reach of these and other topologically identified particle measurements 
(i.e. S(S), and will reach their limit when the momentum of the daughter 

tracks becomes too high to be accurately measured in the TPC. As the momentum 
resolution worsens the invariant mass calculation will become less accurate. As a 
result, the width of the mass peak will broaden. In addition, low pt particles mis- 
measured as high p T particles will start to dominate the less prominent high p T 
signal (feed-down) . The pr scale where the analysis fails has not been extensively 
studied but should depend on the specific particles decay topology. We naively 
expect the Kg identification to fail first, around Pt ~ 15 GeV/c, where the 
high p T signal will be dominated by low p T feed-down. For comparison, the 
7T° identification in the EMC is limited by the detector technology to roughly 
1.5 < p T < 20 GeV/c. 

With detector upgrades and increased data samples, STAR has the potential 
to measure the yield of heavy-flavor mesons and baryons (particularly for D 
mesons), charmonium production (J/V), and direct photon production. Given its 
extensive array of particle identification and event characterization capabilities, 
the STAR detector is particularly well suited for characterizing the matter created 
in heavy ion collisions. 

2.4.1 The STAR Trigger Detectors 

The bunch crossing rate at RHIC is ~ 10 MHz while the read-out rate for the 
STAR TPC is ~ 100 Hz. When the interaction rates approach the bunch crossing 
rates, the STAR trigger must reduce the event rate by five orders of magnitude. 
The STAR trigger needs to reject background, such as beam-gas interactions 
(expected rate ~ 100 Hz), select events that best further our physics goals, and 
issue triggers to the other detectors. Furthermore, the future success of STAR 
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may depend on the ability to trigger on rare events. 

With recent upgrades, the STAR detector system has four fast detectors that 
can be used as trigger detectors: the central trigger barrel (CTB), the zero-degree 
calorimeters (ZDC), a multi-wire counter (MWC), and the barrel electromagnetic 
calorimeter (EMC). In addition, a beam-beam counter (BBC), a forward 7r° de- 
tector (FPD), and an endcap electromagnetic calorimeter (EEMC) will become 
available for triggering |Bie03j . 

The CTB measures the charged particle multiplicity. With 240 scintillator 
slats each covering 7r/30 in <p and 0.5 in rj, the whole CTB covers —1.0 < r\ < 1.0 
and — 7r < < it at a radius of four meters. Its multiplicity resolution is ~ 3% 
for multiplicities > 1000. 

Each RHIC experiment has two ZDC's to monitor beam interactions. The 
ZDC's detect the neutrons freed from the Au ions when a collision occurs (spec- 
tator neutrons). The STAR ZDC's are located ±18.25 m from the nominal in- 
teraction region and subtend an angle 9 < 0.002 radians. Each ZDC consists of 
three modules with a series of tungsten plates and layers of wavelength shifting 
fibers that route Cherenkov light to a photo-multiplier tube. The timing of the 
ZDC signals is also used to locate the longitudinal position of the interaction 
vertex. 



Trigger 


Conditions 


Hadronic Minbias 
Hadronic Central 


[ZDCe > 5 & ZDCw > 5] & CTB > 15 mips 
[ZDCe > 5 & ZDCw > 5] & ZDCsum < 85 
& [Vertex Cut] & CTB > 2000 mips 



Table 2.5: The two trigger settings used in this thesis from the 2001 Au+Au data 
taking. 

During the 2000 and 2001 Au+Au runs the CTB and ZDC were used to study 
minimum-bias, peripheral and central Au+Au collisions. Table I2~5l lists the ZDC 
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and CTB conditions for the two trigger settings used in this analysis; hadronic 
minimum-bias and hadronic central. Figure 12.51 illustrates the selection scheme 
for these triggers in the ZDC verses CTB plane. 
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Figure 2.5: A diagram illustrating the STAR trigger scheme for central (left) and 
minimum-bias (right) triggers. A central event will have a low ZDC count and a high 
CTB count. 



2.4.2 The STAR Time Projection Chamber 

The STAR TPC (Figure I2.6j) surrounds the beam-beam interaction region. The 
inner and outer radii of its drift volume are 50 cm and 100 cm respectively. The 
drift length from the central membrane to either of the ground planes is 209.3 
cm. The central membrane is typically held at 28 kV. A chain of 183 resistors 
and equipotential rings along the inner and outer field cage create a uniform drift 
field from the central membrane to the ground planes where the anode wires and 
pad planes are organized into 12 sectors. 

Figure 12.71 shows a cutaway view of the readout pad planes of an outer sub- 
sector. The first of three wire planes is used as a gating grid. The anode wires 
are located between a shielding wire plane and the cathode pad plane. In the 
open configuration the voltage on the gating grid wires is set so that ions pass 



41 




Figure 2.6: Sectioned view of the STAR TPC showing the inner/outer field cage, the 
central membrane, and inner/outer sectors. 

through freely. When it is closed the field lines terminate on the gating grid wires 
and the electrons and ions cannot pass. When the TPC is not being read-out 
the gating grid is closed and prevents ions from drifting back into the TPC drift 
volume where they can interfere with the uniformity of the drift field. 

The second wire plane shields the TPC drift region from the strong fields 
around the anode wires. As electrons drift past the gating grid and the shield 
plane they accelerate towards the anode wires and initiate a charge amplifying 
cascade. The x-y position of the electron-ion pair left in the TPC by a high 
energy particle is determined by the position of the cathode pads that detect the 
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To I 1 

Figure 2.7: Cutaway of the pad plane of a TPC outer sub-sector. 

cascade. The z position is determined by the time bucket and the drift velocity. 
With 136,608 pad planes and 512 time buckets, the TPC has over 70 million 
three-dimensional pixels. In addition, we use the signal from three adjacent pads 
to better determine the cluster centroid and so, the resolution in the pad row 
direction is significantly smaller than the pad size. The resolution depends on 
the position and orientation of a track relative to the pads, but it is typically 
0.5-1.0 mm. 

Figure 12.81 shows one sector of the TPC pad plane. The inner sub-sector is 
designed to handle the higher track density near the collision vertex. Table 12.61 
lists the dimensions of the inner and outer sub-sectors. Because of the size of the 
front-end electronics, the inner-pad coverage cannot be made continuous. 

In addition to tracking charged particles, the TPC is also able to identify 
particles by their mass. High energy charged particles lose energy as they tra- 
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Figure 2.8: A diagram of a TPC pad plane sector. 





Inner Sub-sector 


Outer Sub-sector 


Pad size (mm) 


2.85 x 11.50 


6.20 x 19.50 


Pad isolation gap (mm) 


0.5 


0.5 


Pad rows 


13 


32 


Number of pads 


1750 


3942 


Anode to pad spacing (mm) 


2.0 


4.0 


Anode voltage (V) 


1170 


1390 


Anode gain 


3770 


1230 



Table 2.6: Geometry of the inner and outer sub-sectors. 



verse the TPC gas. The average energy loss depends on their velocity, not their 
momentum p. At a given p below 0.8 GeV/c pions, kaons and protons suffer sig- 
nificantly different average energy losses. As such, in this region, measurements 
of the energy deposited along a particles trajectory can be used to identify the 
particle. Figure l2~^l shows the energy loss (dE/dx) measured from tracks in the 
STAR TPC where the bands correspond to particles with different masses. 
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Figure 2.9: The energy loss of charged particles traversing the TPC gas. 
2.4.3 STAR TPC Gas System 

The TPC gas system |Kot03j supplies the TPC with either one of two gas 
mixtures — P10 (Ar 90% + CH 4 10%) while the detector is operating or C2H 6 
50% + He 50% for purging the TPC when it is not in use. The TPC gas mixture 
must satisfy multiple requirements. It is the medium where the particles being 
tracked induce ionization, the medium those electron-ion pairs drift through, and 
the medium where the electron multiplication takes place. The convenience and 
safety of the gas is also considered. 

The electron drift velocity in P10 is relatively fast and it peaks and saturates 
at a relatively low electric field (130 V/cm). Operating with a drift field in the 
saturated region minimizes variations in the drift velocity. Some of the important 
characteristics of P10 are listed in table I2~71 
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Drift Characteristics 



Drift Velocity (Maximum) 
Longitudinal Diffusion 
Transverse Diffusion 



5.45 cm//xs 

320 fim/ sqrt(cm) 

185 fim/sqrt(cm) 



at 130 V/cm 
0.5 Tesla Field 
0.5 Tesla Field 



Ionization Characteristics 



227 electrons from a 5.9 keV X-ray (Fe 55 ) 



Charge Created 



Gain (N/Nq) Characteristics 



Inner Sector Gain 


« 3770 


for V anode — 


1170 V 


Outer Sector Gain 


« 1230 


for V anode = 


1390 V 



Table 2.7: Characteristics of the TPC gas mixture, P10. 

The drift velocity and the gas gain are both sensitive to the pressure and 
purity of the gas. The TPC gas pressure varies with atmospheric pressure so 
both of these are monitored. Ionization induced by lasers at fixed locations in 
the TPC are used to measure the drift velocity. The gas gain is monitored by a 
gain chamber and by observing the energy loss of tracks in the TPC. Table EI 
lists characteristics of the TPC gas system. 



System Characteristics 



TPC Volume 
Internal TPC Pressure 
Recirculation Flow 
Oxygen Content 
Water Content 



50000 liters 
2.0 ± 0.03 mbar 
36000 liters/hour 

< 25 ppm 

< 20 ppm 



Table 2.8: Gas system parameters. 



2.4.4 TPC Gas Gain Monitor 

The UCLA nuclear physics group has been responsible for the construction 
and installation of a chamber designed to keep a minute-by-minute record of the 
gain of the STAR TPC gas. At present, we base our corrections for changes in 
the gas gain — thought to be primarily due to variation in the gas pressure — on 
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measurements of the TPC gas pressure. An attempt is also made to correct 
for residual variations in the gain by estimating the ionization energy loss in 
the TPC gas (dE/dx) for tracks which qualify as proton candidates. This step 
requires averaging together data taken over several hours. The final adjusted gain 
is then used to make a better measurement of the dE/dx of tracks as they traverse 
the TPC. As described in Section 12.4.21 this critical measurement is correlated 
with a particles momentum to provide a means of particle identification (PID). 




Figure 2.10: The gain monitor chamber with the wire planes and the shielded elec- 
tronics exposed. The shielding for the wire planes and the Fe 55 source are not shown. 
The aluminum outer housing is folded back and locked into place when the chamber is 
in use so that the high-voltage wires are not exposed. 

With the construction and installation of the new gain monitor chamber com- 
plete, we are able to measure the TPC gas gain directly and as such may be able 
to improve the measurement of dE/dx and PID at STAR. A source of radiation 
with a known energy is used to ionize the TPC gas flowing through the chamber. 
The deposited charge accelerates toward the chamber's high voltage anode wires 
and a charge amplifying cascade ensues. The pulse generated by the cascade 
depends on the energy of the incident radiation and the gas gain. We use an Fe 55 
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source emitting 5.9 keV photons. The signal from the anode wires is amplified 
and conditioned with an Amptek A225 pre-amplifier and shaping amplifier and an 
Amptek A206 voltage amplifier and low-level discriminator. The magnitudes of 
these pulses are analyzed with an Amptek (PMCA600A) multi-channel analyzer 
(MCA). 

We fit the spectrum of pulse heights from the MCA to an exponential function 
for the background noise and two Gaussian functions; one for the 5.9 keV peak 
and another for the secondary photon escape peak at 2.7 keV. The variation in 
the position of the primary peak is used to monitor the relative magnitude of the 
gas gain. We are able to keep the noise level low in the spectrum and have found 
that the resolution of the Gaussian peak is approximately 13%. 

The MCA is read and controlled by a PC (BEATRICE.STAR.BNL.GOV) 
located in the data acquisition (DAQ) room of the STAR hall. No wires can 
connect the 'outside world' to the electronics platform where the MCA is located 
or to the detector where the chamber is located. Instead, the PC controls the 
MCA via optical fibers and a SITECH 2506 fiber-optic modem. The chamber 
is attached to a pipe flange on one of four exhaust manifolds on the face of the 
TPC. Its wire planes extend into the pipe where the P-10 exhaust flows from the 
TPC. The chamber (shown in Figure 12.10)1 is built to replicate the behavior of 
the TPC pad planes. Its geometry — including the diameter of the wire used — 
matches the geometry of the outer sub-sector pads in the TPC. The chamber is 
electrically isolated from the pipe and the wire planes are shielded from stray 
fields by a wire mesh cage surrounding them. The Fe 55 source is mounted inside 
the wire mesh cage several centimeters above the wire planes. When in operation 
the anode wire plane is held at +1390 V. The anode plane is located between a 
grounded pad plane and a wire ground plane. The chambers face is constructed 
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out of non-conducting material so that the bolts used to attach the chamber to 
the exhaust manifold are isolated from the rest of the chamber. This is necessary 
because the exhaust manifold does not share the same ground as the TPC and 
the monitor chamber. 
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Figure 2.11: Left: Gas gain versus anode wire voltage. Right: Pressure dependence 
of the gas gain in P10. 



The chamber was tested at UCLA before installation. The gas gain variation 
with respect to the anode wire voltage is shown in Figure 0. Ill (a). Although 
we do not anticipate that the gain monitor will ever be used with any voltage 
other than +1390 V, we ran this test to compare to other gain monitors and to 
understand how variations in the supply voltage could affect the measured gain. 
Figure 12.111 (b) shows the dependence of the gas gain on the gas pressure. The 
expected relationship of gain to pressure is given by: 



dG dP 
~G ~ ~ A T> 

log(G) w -Alog(P)+B, 



(2.1) 
(2.2) 
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where G and P are the gain and pressure respectively. The coefficient A is 
expected to be 6.7 ±1.5 [BR94 and B is an arbitrary constant. In our calibrations 
we find A = 6.277 ± 0.003 where the error is statistical only. 

The gain monitor recorded data during the entire 2001-2002 data taking. 
The software proved to be very robust and required little or no maintenance or 
intervention from the detector operators or shift crew. Several troubling features 
in the data however, became apparent during the data taking. It was noted early 
in the Run-2 that the gain measurement varied systematically with the magnetic 
field setting (shown in Figure |2\12|) . Possible causes include variations in the 
response of the electronics with the magnetic field setting, a change in the charge 
amplification caused by the orientation of the wire planes relative to the magnetic 
field, or gas leaks exacerbated by the magnetic field. For the 2003 d+Au collisions 
the gain monitor was realigned to match the orientation of the TPC pad planes. 

We also find that the relative magnitude of the gas gain decreased with time; 
as seen in Figure l2".13l fb). This time dependence was observed in other calibration 
data and has been accounted for in the dE/dx calibrations. The importance of the 
gain monitor chamber however, can be seen in Figure 0.131 (a). Scatter is seen in 
the plot of log(G) versus log(P). This may indicate that there are still variations 
in the gain not taken into account with the current calibration method — a method 
that is insensitive to gain variations on a short time scale. 

Further study of the gain monitor and the gain monitor data is needed before 
we use it for calibrations. The electronics were returned to UCLA where they 
were tested in a 0.3 Tesla magnetic field for variations in pulse height with field 
direction. No effect was observed. We have also built and installed an adapter 
that rotates the gain monitor so that the wire planes are perpendicular to the 
magnetic field. In the 2001-2002 data taking the gain monitor data was recorded 
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Figure 2.12: The variation of the gain with the magnetic field — after accounting for 
the pressure dependence — is shown for four different time periods. The dependence 
is not strong in the upper-left plot (earlier in the Au+Au Run-2) but subsequent 
measurements show a strong field dependence. In the inset of panel (d) the vertical 
axis is adjusted to emphasize the field dependence. 

in the online database, but was not propagated to the off-line database where it 
can be easily used for calibrations. This will be changed for future data taking 
and we anticipate that the gain monitor will be used for calibrations as well as 
diagnostics. 
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Figure 2.13: Left: Scatter plot of log(G) versus log(P). The contours represent the 
expected behavior for log(G) = — 6.277 log(P) + B. Right: Variation in gain with time: 
the mean value of the "constant" B is shown after the expected pressure dependence 
has been parameterized. 
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CHAPTER 3 



Analysis Methods 

The technique for finding K s , A, or A candidates with the STAR detector and 
calculating their p T and y distributions is well established |Adl02e| lLon02j 1 . Up 
to now, measurements of identified particle t>2 have relied on pure particle iden- 
tification (90% purity 2 ) via dE/dx measurements in the TPC gas |Adl01j . We've 
adapted the Vi analysis method to calculate v% for particles identified only on 
a statistical basis. With this method it is possible to calculate t>2 for identified 
particles independent of the particle sample's signal-to-background ratio. 

In this chapter, we describe the selection criteria for events, tracks and the K s , 
A, or A candidates. Details of the Rcp and spectra measurements — including an 
analysis of systematic errors — are given. Finally, the analysis methods for mea- 
suring t> 2 of K s , A, and A are presented along with a discussion of the systematic 
errors associated with the Vi analysis. 

3.1 Event and Track Selection 

To date, for the STAR experiment, the number of events that are useful for 

lr The author thanks H. Long for his assistance in making the measurements presented in 
this thesis. His work with weak-decay-vertex finding has become a cornerstone of the STAR 
collaborations scientific program. 

2 Purity is defined as the raw yield of the particle at a given dE/dx value, divided by the 
sum of all other particle yields with the same dE/dx. 
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Data set 


Minimum-bias 


Central 




Recorded Used 


Recorded Used 


Kun-1 (^/s NN 
Run-2 Qs NN 


= 130 GeV) 
= 200 GeV) 


0.8 x 10 6 0.20 x 10 6 
4.0 x 10 6 1.6 x 10 6 


0.8 x 10 6 0.18 x 10 6 
3.5 x 10 6 1.5 x 10 6 



Table 3.1: The number of Au+Au events used in our analysis. 



our analysis has been ~ 40% of the total recorded. Table 13.11 lists the number 
of events recorded and used for Run-1 and Run-2. Events for which no primary 
vertex is found are discarded. For Run-1, events with z- vertex further than 75 cm 
from the TPC center were discarded. For Run-2, improvements in the accelerator 
allowed STAR to select only collisions within 25 cm of the TPC center. Still 
more events are discarded to remove trigger biases. A large sample of minimum- 
bias data was taken with a tight ^-vertex cut applied in the level zero trigger 
that biased the sample. These events require more careful analysis and are not 
included in this analysis. 




Figure 3.1: The TPC charged particle multiplicity distribution for ^s NN = 200 GeV. 
N c h is the number of primary tracks in \r/\ < 0.5. The fractions of the total cross-section 
used for the analysis of v<i are shown. 
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The multiplicity as measured by the TPC — not the CTB — is used to define 
STAR's centrality intervals. The TPC reference multiplicity for Run-2 is the 
total number of primary tracks in the TPC with 10 or more fit points, having 
\rj\ < 0.5, and a distance of closest approach (DCA) to the primary vertex less 
than 3 cm. A primary track is defined by a helix fit to the TPC points and to 
the primary vertex; the global track fits do not include the primary vertex. For 
Run-1 primary tracks within \rj\ < 0.75 were used to define the multiplicity. 



Event plane K° s and A (A) 
Track set Primary Global 
DCA to primary vertex (cm) < 3 na 
Number of hits > 15 > 15 

Number of hits/possible hits > 0.52 na 

\q\ < 1.5 na 



Momentum (GeV/c) 0.1 < p T < 2.0 p T > 0.05 

Table 3.2: Track selection criteria for Au+Au collisions at y/s NN = 200 GeV. 

Table 13.21 list the selection criteria for tracks used in the analysis of y/s NN = 
200 GeV data. For the K s , A, or A reconstruction when the dE/dx of a track 
can be used to identify the particle type, an additional dE/dx cut is made. For 
0.2 < pt < 2.0 GeV/c, pion candidates are required to have a dE/dx value within 
6 < cr,,- < 5 and proton candidates are required to have a dE/dx value within 
2.85 < ij p < 10. These cuts are very loose and only act to exclude tracks which 
are obviously not of the correct type. The most effective selection criteria in the 
identification of K s , A, or A particles are the decay topology cuts. 

3.2 Decay Vertex Topology: Yield Measurements 

We identify the K s , A and A candidates from the charged daughter tracks 
produced in the weak decays: K s — > 7r + 7r~, A — > pn~ and A — > pn + . To 
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Figure 3.2: A sketch of the geometry of a neutral particle decay in the TPC. Case 
"a" shows charged daughter tracks curving towards each other after they are created 
in the decay while in case "b" they are curving away from each other. 
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select the K s , A, or A candidates, we calculate the distance of closest approach 
(DCA) between all combinations of selected global tracks within an event. We 
define the four-momenta of the daughter particles by assuming they originated 
from the points on the two helices where the DCA occurs, and by choosing 
a mass hypothesis appropriate for the weak-decay channel. We use the four- 
momentum of the two daughter particles to calculate the invariant mass and 
kinematic properties of the candidate. 



Candidate px (GeV/c) 


< 


1.6 


1.6-3.0 


> 


3.0 


Daughter-daughter DCA 


< 


0.80 


< 0.80 


< 


0.80 


Daughter- V prim DCA 


> 


3.00 


> 3.00 


> 


0.50 


Decay Length 


4.0 


-25.0 


4.0-40.0 


5.0 


-60.0 


V — Vprim DCA 


< 


0.80 


< 0.80 


< 


0.80 



Table 3.3: Selection criteria for the final analysis of the K s Rcp- Units are centime- 
ters except where indicated. 

Further selection criteria (i.e. cuts) are applied to the orientation of the two 
tracks — with respect to each other and with respect to the primary vertex — to 
increase the probability that the track combination is associated with a real decay. 
Figure 13.21 illustrates the geometry of a neutral-particle decay vertex (V°) . For 
Kg, A or A decays in a magnetic field, two equally probable cases occur: the 
daughter tracks curve towards each other or the daughter tracks curve away from 
each other. The geometric variables used to select Kg, A or A decays are shown 
in the figure. Table l3~3*l shows the K s selection criteria used for the spectra and 
Rcp analysis. We choose the vertex geometry cuts to minimize the statistical 
and systematic uncertainty in the measured Rcp- 

3.2.1 Invariant Mass Distributions 

The Kg, A, or A particles are not identified on a particle-by-particle basis but 
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their uncorrected yields are extracted from the peak at their known masses in the 
invariant mass distributions. The yield is estimated by fitting a smooth function 
to the combinatorial background outside the peak region. We determined that 
the background is dominated by combinatorial counts by rotating all positive 
tracks 180 degrees in the transverse plane and reconstructing the K s and A (A) 
decay vertices. This procedure destroys all real vertices within our acceptance 
so that we can describe the combinatorial contribution to the invariant mass 
distributions. 




Figure 3.3: Invariant mass distributions for ir + ir from 0-5% central collisions (top) 
and 40-60% central collisions (bottom) in three px intervals. 

The observed masses, 496 ±8 MeV/c 2 for 7r + 7r~ and 1116 ±4 MeV/c 2 forp7r, 



58 



are roughly consistent with accepted values [GGOOJ and the widths are determined 
by the momentum resolution of the detector. For px < 1.5 GeV/c, however, the 
Kg peak is shifted to a lower mass. At px = 0.2 GeV/c the peak is shifted by 
the greatest amount, 10 MeV/c. This shift is, for the most part, replicated by 
simulations and is attributed to energy loss suffered by the daughter particles in 
the detector material. Figure 13.31 shows invariant mass distributions for the K s 
Rcp analysis. When the same selection criteria are used for all px and centrality, 
the combinatorial background is larger for lower px and for more central events. 
At higher px, where the size of our data sample is limited, we place less stringent 
requirements on the candidates. As a result, the combinatorial background is 
quite large in central events for px > 3.0. 

3.2.2 Detector, Tracking and Reconstruction Efficiency 

Simulations are used to calculate the efficiency of the detector and the tracking 
software |Lpn02j. The TPC response to Monte-Carlo generated Kg, A, or A 



decays is simulated. The simulated clusters (the pixel level TPC response) are 
then embedded into real events and these events are passed into the Kg, A, or 
A reconstruction chain. Reconstructed candidates are then associated with the 
embedded particles so that the efficiency of the detector and the reconstruction 
chain can be estimated. 

Figure EOl shows the efficiency correction factor for Kg in the three centrality 
intervals used in the calculation of Rcp- The number of simulated K s particles 
we embed at a given px is varied to approximate the true px dependence of 
the yield. The slope of the px spectrum is characterized by the inverse slope 
parameter T of an exponential fit. By matching the slopes in the simulations 
to the real slopes, the bin-sharing for the simulated particles replicates the bin- 
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Figure 3.4: The -Kg efficiency for the selection criteria listed in Table 13.31 The 
discontinuities in the efficiency at pt = 1.6 and 3.0 GeV/c reflect changes in the 
selection criteria (see Table l3~3"]) . The line histograms show the relative magnitude of 
the statistical error for the efficiency 

sharing for real particles. In this way, the efficiency correction also acts as a feed- 
down correction. With a slope on the embedded particle yield, it is necessary to 
embed within limited pt ranges to generate high statistics at large pr- To match 
the feed-down in real data, simulated data cannot be used near the edge of the 
embedded px range. Table l3~4l lists the slope parameters, the pt ranges, and the 
number of events used for the embedding. 



Embedded range (GeV/c) T (MeV/c 2 ) px coverage (GeV/c) Events 



0.0-1.4 300 0.0-1.2 39 k 

1.0-2.2 315 1.2-2.0 45 k 

1.8-2.8 330 2.0-2.6 49 k 

2.4-3.4 330 2.6-3.4 36 k 

3.0-5.2 450 3.4-5.0 15 k 

4.6-10.0 500 5.0-8.0 19 k 



Table 3.4: Embedded data for K s analysis. 



3.2.3 Systematic Uncertainties 

Systematic errors for the spectra and Rep are introduced from uncertainties 
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in the detection efficiency, the reconstruction efficiency, the background subtrac- 
tions and from mis-measurements of the candidates px- The momentum resolu- 
tion Spt/pt is estimated from simulations. The other systematic uncertainties 
are studied by varying the Kg, A, and A selection criteria. Changing the se- 
lection criteria varies the number of background counts and tests how well the 
distributions in the simulated data match the real data. When all the relative 
distributions are accurately simulated, changing the cuts will not change the 
efficiency corrected particle yields. 



T 

, i , i , 

2 4 6 

^-Vprim DCA CUt ( Cm ) 

Figure 3.5: The K s yield for different it^ 1 — V vr im DCA selection criteria, relative to 
a 3.0 cm cut. 

We found that the K s yield depends strongly on the ir^ — V pr i m DCA cut. 
This variation is most likely caused by either space-charge or ExB distortions 
which arise as ionization drifts toward the TPC pad planes. Figure 13.51 shows 
the measured K s yield for several n ± — V pr i m DCA values. The yields are scaled 
by the yield from the tt ± — V pr i m DCA > 3.0 cm analysis. No px dependence 
is seen in the variation of the scaled yield, so the px integrated yield is used in 
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Figure 1331 For 3.0 cm and above the yield is independent of the n ± — V pr i m DCA, 
so we used a 3 cm cut for our final analysis. 
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Figure 3.6: Kg yield for forward and reversed field settings. 

Figure EOH shows the Kg spectra for forward and reversed field settings. Early 
in the analysis an approximately 8% systematic deviation was seen between the 
yields from events with the two field settings. By selecting only Kg candidates 
with a decay length from 4-20 cm the deviation is removed. Final estimates 
for the systematic errors on Rcp are dominated by variation in the yields with 
different selection criteria. Table I3.5I lists systematic uncertainties for Kg and 
A + A R C p. 
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A ± A 


p T (GeV/c) 


1.0 2.5 4.0 


1.0 2.5 4.0 


Rcp (bg) 
Rcp (eff) 


±0.04 ±0.02 ±0.08 
±0.10 ±0.10 ±0.10 


±0.02 ±0.04 ±0.06 
±0.10 ±0.10 ±0.10 


5p T /p T 


±0.012 ±0.027 ±0.030 


±0.016 ±0.027 ±0.037 



Table 3.5: The systematic errors from background (bg), and the efficiency calculation 
(eff) are listed for Rcp (0-5%/40-60%) at three pt values along with the pr reso- 
lution (Spt/pt)- The values listed are relative errors and do not include the overall 
normalization uncertainty in Rcp from the calculation of Nbi n - 

3.3 Reconstructing the Reaction Plane 

The real reaction plane is not known, but the event plane, an experimental 
estimator of the true reaction plane, can be calculated from the azimuthal dis- 
tribution of primary tracks |PV98j . The select ion criteria for the primary tracks 
used to calculate the event plane are given in Table l3~2l We require the ratio of 
the number of space points to the expected maximum number of space points 
for each track to be greater than 0.52 to prevent split tracks from being counted 
twice. For the analysis using y/s NN = 130 GeV data the events are required to 
have a primary vertex within 75 cm longitudinally of the TPC center (z- vertex). 
During the Au±Au running with ■ s /s NN = 200 GeV the ^-vertex distribution was 
narrower so those events are required to have |z-vertex| < 25 cm. These cuts do 
not bias our analysis. 

For experiments at RHIC energies — unlike those at lower energies — the reac- 
tion plane is assumed to be transverse to the beam axis. As such, since it is not 
necessary to rotate the flow coordinate system in the polar direction, only the 
transverse direction is considered. The error introduced by this assumption will 
go as the square of the polar flow angle, 0fi ow ~ (Px) /Pbeam <C 1, and is negligible. 

With a perfect detector the azimuthal distribution of the event plane would 



63 



be isotropic. In a realistic environment, however, limitations to a detectors ac- 
ceptance lead to a bias in the estimation of the reaction plane. Acceptance 
corrections are introduced to account for both the limited coverage and the im- 
perfect efficiency of a real detector. With the STAR detector, detector biases are 
removed from the event plane distribution by applying weights to the tracks used 
in its calculation. The 0-weights are generated by inverting the normalized </> 
distribution for tracks from many events. When other weights are included in the 
event plane calculation (e.g. pt weights) they are also applied to the 0- weights. 
In this way, after the 0-weights are applied the azimuthal distribution for tracks 
is — by construction — isotropic. 

The acceptance of the STAR detector system depends on the kinematic vari- 
ables of a particle (77, <fi and pr), the longitudinal position of the collision vertex 
in the TPC, the multiplicity of the event and the hour-by-hour state of the de- 
tector. During Run-2 electronics failures resulted in the temporary removal of 
read-out (RDO) boards from the data acquisition chain. The masking and un- 
masking of RDO boards changed the detectors acceptance with time. To ensure 
that all these variations are accounted for properly, the 0-weights are calculated 
separately for positive or negative rj, for positive or negative z-vertex position, 
for magnetic field polarity, for nine different centralities and for four different 
detector states. 

Figure 13.71 shows examples of the 0-weights. The 12-fold periodic structure 
is caused by the change in the detectors acceptance near the sector boundaries. 
The 24-fold periodic structure arises because a track that starts near the edge 
of a sector and curves toward the middle of the sector has a high probability 
of being reconstructed. As a consequence, the efficiency for detecting a positive 
particle is enhanced near the one side of a sector while the efficiency for detecting 
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Figure 3.7: An example of the weights used to correct the reaction plane calculation 
for acceptance is shown. These weights are for the centrality bin corresponding to 
30-40% of the Au + Au collision cross-section and forward field polarity with the RDO 
board 3 in sector 9 and RDO board 4 in sector 21 where masked out. 

a negative particle is enhanced on the other side: two maxima are seen for each 
sector. 

When the tracks in the event plane calculation are re-weighted properly the 
distribution of the azimuthal angle of the event plane is isotropic. Figure 13.81 (a) 
shows the event plane's azimuthal distribution fit to a constant. When the event 
plane distribution is flat, the acceptance for the particles being compared to the 
event plane will not bias the measurement of v-i'- as long as either the distribution 
of the event plane or the distribution of tracks is isotropic, acceptance effects 
introduce no bias. Poor acceptance will, however, negatively impact the reaction 
plane resolution. The decrease in the reaction plane resolution caused by the 
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imperfect acceptance can be accounted for by applying a resolution correction 
factor to the measured v^- 
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Figure 3.8: Left: Azimuthal distribution of the event plane fit to a constant. The 
event plane distribution should be 'flat' (azimuthally isotropic). Right: Distribution 
of the difference between the event-plane angles for sub-events of randomly assigned 
tracks. 



A correction to the observed Vi is introduced to account for the uncertainty in 
the determination of the reaction plane (i.e. the reaction plane resolution). 
The event plane (calculated from Equation 13.4)) is used to estimate the reaction 
plane: the accuracy of the estimation depends on the number of tracks used 
and the magnitude of the true v-i signal. With an infinite number of tracks, a 
finite v 2 , and a perfect detector the event plane could be a perfectly accurate 
estimator of the true reaction plane. With a limited number of tracks detected 
in a real detector, however, we cannot assume that the true reaction plane has 
been accurately estimated. Poor resolution leads to a decrease in the calculated 
value of the flow parameters because, the correlation of the particles with the 
reaction plane is partially lost. The correction factor necessary to compensate 
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for the resolution is found as follows: 



v * = ( e ai(*-*«,)) jV ok = ( e 2i(<k-* BV )) (3.1) 
\ = (exp 2 ^--*^)), (3.2) 

where vf, v% bs , ^%p, and are the real V2, the observed t>2, the real reaction- 
plane angle, and the reconstructed event-plane angle. From Equation 13.21 the 
proper correction factor is found to be (cos2( 1 3/ ep — ^rp))- This quantity can 
be calculated by reconstructing event planes from two random sub-sets of tracks 
within the same event (sub-events). The difference between the two sub-event- 
plane angles ($ EP — ^ep) * s shown in Figure (b). The resolution measured 
from two sub-events with equal numbers of tracks is given by Equation 13.31 



(cos [2 (*i P - ¥«,)] ) = ^(cos [2 (*j. p - ¥*,)]). (3.3) 

The second harmonic sub-event-plane angles are calculated from Equation 13.41 

/ t x Y] s Wi sin (26;) 

tan {2^ EP ) = ^ (3.4) 

2^ wi cos (2(pi) 

The tOj's in Equation 13 .41 are weights used to maximize the resolution. In our case 
we use the particle's transverse momentum and the 0-weights as the weighting 
factor. An interpolation formula is then used with an iterative routine to calculate 
its roots and find the correction factor for the full event plane {cos2(ty ep~ ^%p))- 
We estimate the errors in the correction factor by varying the input by a small 
amount and calculating the change in the result. 

The resolution correction factor is shown in Figure 13.91 The resolution de- 
pends on the number of tracks used and the magnitude of the event asymmetry. 
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Figure 3.9: Left: Resolution correction factor for the y/s~^ = 200 GeV analysis. 
Right: Resolution correction factor for the y/s NN = 130 GeV analysis. The increase 
in the resolution from Run-1 to Run-2 reflects the implementation of a new method 
(discussed in the text) for calculating the event plane. The resolution for the Run-1 
A analysis was lower than the K s analysis, only reaching a maximum of 0.58 ± 0.007. 
With the new event plane calculation the resolution for A and K s v 2 are the same and 
reach a maximum of 0.813 ± 0.003. 

For the most peripheral events the small number of tracks available reduces the 
resolution while for the most central events the symmetry of the collision overlap 
region degrades it. As a result, the resolution is greatest at a centrality corre- 
sponding to roughly 20-30% of the collision cross-section. 

Since t> 2 is calculated from the distribution of (0j — if particle % is in- 

cluded in the evaluation of ^rp, auto-correlations are introduced. For the 
measurement of Kg v 2 in Run-1, these unwanted correlations are eliminated by 
calculating the event plane from charged particle tracks with a DCA to the pri- 
mary vertex less than 1 cm while we use only tracks with a DCA greater than 
1 cm to reconstruct the K s . For A we exclude from the event plane calculation 
all tracks identified as proton candidates by their energy loss in the TPC. This 
method has one important disadvantage: the reduction in the number of tracks 
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used degrades the resolution of the reaction plane. This in turn leads to larger 
errors on the final t>2 measurement. 

To maximize the number of tracks available for our event plane calculation, 
we implemented a new method to remove auto-correlations for neutral vertex 
particles. Rather than dividing tracks into two sets — one to be used for re- 
constructing neutral vertices and the other to be used for calculating the event 
plane — we calculate a separate event-plane angle for every K s , A, or A candidate. 
Auto-correlations are removed by excluding only the two tracks associated with a 
specific K s , A, or A from the event plane calculation. With a larger track sample 
the reaction plane resolution increases and the statistical and systematic uncer- 
tainty on v 2 decreases. Figure 13.91 shows the reaction plane resolution correction 
factor from the method used for the y/s NN = 130 GeV data and the method used 
for the y/s NN = 200 GeV data. 

3.4 Calculating the v 2 of K Q S and A + A 

We use the measured yield within multiple intervals of (fa — ^ l RP ) to calculate 
t>2 = (cos [2(fa — ^ijp)]}, where fa is the azimuth angle of the momentum vector 
of particle i and ^ RP is the reaction-plane angle for the event that particle i 
was observed in. To remove autocorrelations, the contributions from the decay 
daughter tracks associated with particle i are subtracted from the right hand 
side of Equation 13.41 The candidates are categorized by invariant mass, px and 
(fa — ^ijp). For each px interval, twenty ((f) — ^ RP ) intervals from to 2n are 
created (all events are combined). The yield in interval j (drij) is calculated by 
fitting a smooth function to the mass region outside the candidate mass window 
and integrating the number of counts above this background. To minimize the 
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systematic errors associated with fitting the background, the same background 
shape is used for every j (0 — ty RP ) interval: only the relative amplitude of the 
background is allowed to change. Once the yields are known, we calculate ^(pt) 
using 




Figure 3.10: Distribution of K s candidates in the invariant mass versus (0 — ty RP ) 
plane. 

Figure 13.101 shows the px inclusive distribution of the Kg candidates in the 
invariant mass versus (0 — ^rp) plane. The V2 signal — an enhancement near 
zero, 7r and 2ir radians — is clearly visible. We find our method of calculating V2 
is insensitive to the background contamination, so, to maximize our statistical 
sample, we use relatively loose selection criteria for the identification of K s , 
A, and A candidates. Table 13.61 lists the criteria used for the final analysis. 
Figure 13.111 shows the same distribution for A and A candidates. 
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Figure 3.11: Distribution of A and A candidates in the invariant mass versus 
((f) - q> RP ) plane. 

3.4.1 Systematic Uncertainties and Correlations Unrelated to the Re- 
action Plane 

Sources of systematic error in the calculation of v-i are correlations unrelated 
to the reaction plane (non-flow effects), estimation of the yield from the invariant 
mass distributions, the pr resolution (Spr/pr), and biases introduced by the cuts 
used in the analysis. Table 13.71 lists the dominant systematic errors for Kg and 
A + A t> 2 . The systematic error in t> 2 associated with the yield extraction is found 
to be small and the non-flow systematic error is dominant. 

The magnitude of charged particle v-i absent of non-flow effects has been 
estimated using a four-particle cumulant analysis [Adl02c : a method thought to 
be less sensitive to non-flow correlations. Figure 13.121 (a) shows the ratio of v 2 
from a four particle cumulant analysis and a two particle cumulant analysis. The 
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(A) 


Pt V / CJ 


^ n a of. on \ n 


/ o n \ 9 n 


7T + (p)-7T~ DCA 

7T+(p)-y pri m DCA 
7T —Vprim DCA 

Decay Length 

V° — Vprim DCA 


< 0.70 < 0.70 < 0.80 

> 1.50 > 1.50 >0.35 

> 1.50 > 1.50 >0.35 

> 3.50 > 5.50 > 7.00 

< 0.70 < 0.70 < 0.80 


< 0.70 < 0.70 

> 0.50 > 0.25 

> 1.00 > 1.00 

> 4.50 > 4.50 

< 0.70 < 0.70 



Table 3.6: The vertex and daughter track selection criteria for Kg, and A. For 
A criteria are the same as A with the p(tt~) exchanged with p(tt + ). All units are 
centimeters unless indicated otherwise. 







A + A 


p T (GeV/c) 


1.0 2.5 4.0 


1.0 2.5 4.0 


V2 (bg) 
v 2 (n-f) 


+0.000 +0.001 +0.003 
+0.001 -0.007 -0.018 
+0.00 +0.00 +0.00 
-0.01 -0.04 -0.03 


+0.001 +0.005 +0.005 
-0.007 -0.001 -0.001 
+0.00 +0.00 +0.00 
-0.01 -0.04 -0.04 


5p T /p T 


±0.012 ±0.027 ±0.030 


±0.016 ±0.027 ±0.037 



Table 3.7: The systematic errors from background (bg) and non-flow effects (n-f) are 
listed for v 2 (0-80%) at three px values along with the px resolution (5pt/pt)- The 
values listed are absolute errors. 

two particle cumulant v 2 result has been shown to be similar to — but slightly 
larger than — the v 2 from a reaction plane analysis [Adl02c . Figure IB*. 121 suggests 
that for minimum-bias collisions, non-flow correlations may account for 10-20% 
of the charged particle v 2 . 

The four-particle cumulant method can be adapted to study v 2 for A + A and 
Kg absent of non-flow contributions but, to be decisive, will require a large data 
sample. Nuclear modification of jet production and fragmentation could lead to 
a particle-type dependence in the relative fraction of the non-flow contribution 
to v 2 . At px > 3 GeV, jet production is thought to be a likely source of non- 
flow correlations. The effect of standard jet fragmentation on v 2 was examined 
using superimposed p+p collisions generated with PYTHIA SjoOl . Within the 
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measured pt region, no significant difference is seen between A + A and K s non- 
flow from this source. As such, in this analysis, we assume a similar magnitude 
for the non-flow contribution to the v 2 of A + A and K s . 
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Figure 3.12: Left: Ratio of V2 calculated by the fourth-order cumulant method and 
the reaction plane method |Adl02c| . Right: Integrated V2 as a function of centrality 
for each method |Adl02c| . 

Figure 13.121 (b) shows the integrated V2 versus centrality (N c h/N max ) from 
Reference jAd!02c| . where v% calculated using a four particle cumulant analysis 
is compared to v 2 from a reaction plane analysis. Table 13.81 lists the values of 



N c h and the collision cross-sections corresponding to the x-axis in Figure 13.121 
(b) |Adl 02cj. The difference between the V2 calculated from these methods is 
used to estimate the centrality dependence of non-flow effects. The cumulant 
analysis indicates that non-flow effects are largest in the most central and most 
peripheral events. Some or all of the difference between these methods could, 
however, arise from event-by-event fluctuations in the initial source shape: these 
fluctuations would reduce the value of t>2 calculated from the cumulant analysis. 
As such, the cumulant analysis and the reaction plane analysis are often taken 
as estimates of the upper and lower limits on the true v%. 
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■N c h/N max 


0.849 


0.708 


0.590 


0.472 


0.363 


0.258 


0.160 


0.060 


X-section (%) 


0-5 


5-10 


10-16 


16-24 


24-31 


31-41 


41-53 


53-77 



Table 3.8: Percent of the collision cross-section corresponding to N c h/N max . The 
value of N max is approximately 878. 
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CHAPTER 4 



Results 



As shown in Figure 1231 K s , A and A particles were identified during Run-2 
across a broader px range than any other particle. In this chapter we present the 
measurement of v-i for Kg and A + A at mid-rapidity from Au+Au collisions at 
y/s NN = 130 1 and 200 GeV. The pr spectra are shown for 0-5%, 40-60%, and 
60-80% centrality intervals. The centrality dependence is studied via the nuclear 
modification factor Rcp which is derived from the spectra. 

4.1 Elliptic Flow 

Elliptic flow at mid-rapidity as a function of transverse momentum for minimum- 
bias (a), 0-11% central and 11-45% central (b) Au+Au collisions at y/s NN = 
130 GeV is shown in Figure |4~T1 The vi of K s and the v-i of A + A both increase 
monotonically with in the 11-45% centrality interval. Throughout the mea- 
sured pt range, the v% for both particles is larger in more peripheral collisions 
than in the central collisions. A similar dependence was observed for charged 
particles in Au + Au collisions at the same RHIC energy jAckOll. Also shown 
in the figure is ^(pt) for charged hadrons Adl03aJ and hydrodynamic model 

1 The v% measurements at ^/s NN — 130 GeV were made in collaboration with J. Fu and were 
published in References |Sor02al IAdl02al ISor02bl lFuB2] . 
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Figure 4.1: Elliptic flow for if<j and A + A particles at mid-rapidity in Au+Au 
collisions at y/s NN = 130 GeV. 

calculations. Within statistical uncertainty, the minimum-bias K s results are in 
agreement with the v-i of charged kaons (not shown) [AdlOlJ. We observe that v-i 
for both strange particles increases as a function of pt up to about 1.5 GeV/c, 
similar to the hydrodynamic model prediction. For p? > 2 GeV/c however, 
the values of t>2 seem to be saturated. It has been suggested that the shape and 
height of t>2 above 2-3 GeV/c in a pQCD model is related to energy loss in an 
early, high-parton-density, stage of the evolution [GVW 01 . 

These are the first measurements of v 2 for K s and A + A at RHIC energy 
and the first measurement of V2 for any identified particle above px ~ 1.0 GeV/c. 
The statistical sample for Run-1, however, is limited and several important de- 
tails remain to be studied — how well do the hydrodynamic models reproduce the 
mass dependence for K s and A + A, does V2 for all particles saturate with the 
same magnitude at the same p?, and does the relative strangeness content of 
the particle affect its elliptic flow? For Run-2, improvements were made to the 
analysis technique and a larger data set became available. 
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eg 




Transverse Momentum p T (GeV/c) 

Figure 4.2: The minimum-bias (0-80% of the collision cross-section) V2(pt) for Kg, 
A + A and at mid-rapidity in Au+Au collisions at ^s NN = 200 GeV. The error bars 
shown are statistical only. Hydro dynamical calculations of V2 for pions, kaons, protons 
and lambdas are also plotted HKH01 . 
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Figure 14.21 shows minimum-bias v 2 at y/s NN = 200 GeV for K s , A + A and 
charged hadrons h + + h~ . The analysis of the charged hadron v 2 is described 
in Reference |Adl03aj . The curves in the figure represent hydrodynamic model 
calculations of v 2 for pions, kaons, protons, and lambdas [HKH01J. At low px, 
the model calculations are in good agreement with the mass and pt dependence 
of v 2 . At intermediate pt however, we find v 2 a+a > v %k in contradiction to 
hydrodynamical calculations: where at a given p?<, heavier particles have smaller 
v 2 values. The p^-scale where v 2 deviates from the hydrodynamical predic- 
tion is ~ 2.5 GeV/c for A + A and ~ 1 GeV/c for K s . Our measurement at 
y/s NN = 200 GeV establishes the particle-type dependence of the v 2 saturation 
at intermediate pr (1.5 < pr < 4.0 GeV/c). 
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Figure 4.3: The v 2 of K s and A + A at ^/s~^ = 200 GeV as a function of pr 
for 30-70%, 5-30% and 0-5% of the collision cross-section. The error bars represent 
statistical errors only. The non-flow systematic errors for the 30-70%, 5-30% and 0-5% 
centralities are -25%, -20% and -80% respectively. 

Figure 14. HI shows v 2 of K s and A + A at mid-rapidity for Au+Au collisions at 
^/s NN = 200 GeV as a function of pr for three centrality intervals: 30-70%, 5- 
30%, and 0-5% of the geometrical cross-section. The px dependence of v 2 from all 
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three centrality bins has a similar trend: a monotonic rise with px at low pr and a 
saturation at intermediate pt- The saturation of vi in the intermediate pt region 
from minimum-bias trigger data in Figure 14.21 is not due to the superposition of 
drastically different px dependencies for various centrality bins. The values of V2 
at saturation show a particle-type and centrality dependence. 



4.2 Spectra 

In Figure 14.41 we show the spectra for K s and A + A [Lon03j at mid-rapidity 
in Au+Au collisions at v /s^7 = 200 GeV. At p T ~ 1.5 and 4.5 GeV/c, the K s 
and A + A yields coincide. In the intermediate pt region, a greater number of 
A + A particles are produced than Kg particles. The K s spectra show a clear 
hardening at higher px while the A + A spectra appears to remain soft within 
most of the measured pt range. As a result, the number of Kg particles produced 
becomes larger than the number of A + A particles again for p? > 4.5 GeV/c. For 
peripheral collisions the separation at intermediate px between the K s and A + A 
yields appears much smaller than in central collisions. The nuclear modification 
factor R C p is a useful measure for studying the relative centrality dependencies. 



4.3 Nuclear Modification Rqp 



Figure l4~5l shows Rcp for K s , and A + A using the 5% most central collisions, 
normalized by peripheral collisions (40-60% and 60-80%). For reference, the 
charged kaon and charged hadron Rcp are also shown. For charged hadrons, these 
40-60% and 60-80% bins have been shown to approximately follow binary collision 
scaling (relative to p+p collisions) without medium modification [Ada03bJ. The 
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Figure 4.4: Spectra for K® and A + A |Lon03| at mid-rapidity from central and 
peripheral Au+Au collisions at y/s NN = 200 GeV. 
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Transverse Momentum p T (GeV/c) 

Figure 4.5: Nuclear modification of K$ and A + A production at mid-rapidity for 
Au+Au collisions at ^/s NN = 200 GeV. For comparison, the Rcp for charged kaons is 
also shown where the kaons have been identified from single-prong decays (kinks) as 
discussed in Reference IAdl02dj . 
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bands in Figure I4~B1 represent the expected values of Rcp for binary and partici- 
pant (N part ) scaling including systematic variations from the calculation [Acla03bJ. 
Table HH] gives the values for Nbm and N part 



Cross-section (%) 


0-5 


40-60 


60-80 


(N bin ) 

(N part ) 


990+$ 
352^ 


91.8 + | 
61.0+i° 


20.0+^ 
19.8+J 



Table 4.1: Monte-Carlo Glauber model calculations of the number of participat- 
ing nucleons N part and the number of binary nucleon-nucleon collisions Nbin for three 
centrality intervals |Ada03b| . 

The kaon and (anti-)lambda yields are suppressed by different magnitudes and 
the pr-scales associated with the onset of the high p? suppression are different. 
Rcp for kaons increases with p? from the participant scaling limit, reaches a 
maximum of approximately 0.6 at pt ~ 1-6 GeV/c and then decreases with px- 
The Rcp for A + A, however, rises to a maximum of approximately 0.9 at px ~ 
2.0 GeV/c, remains near that value up to pt ~ 3.5 GeV/ c and then decreases with 
p T . For most of the intermediate px region, A + A Rcp within errors coincides 
with binary collision scaling, while the kaon Rcp is significantly below unity. For 
both species, the pr where Rcp begins to decrease approximately coincides with 
the pt where t>2 in Figure l4~2l saturates. At highp^ (pt > 5.0 GeV/c), Rcp values 
for Kg and A + A are approaching the value of the charged hadron Rcp- The 
apparent disappearance of the particle-type dependence of Rcp may signify that 
single parton fragmentation dominates the features of Rcp above pt ~ 5 GeV/c. 
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CHAPTER 5 



Discussion 



Event-by-event azimuthal anisotropy in particle production is thought to 
probe the early stages of relativistic heavy- ion collisions |Q1192| . For pt < 
1.0 GeV/c, hydro dynamic models [HKHOlj describe v 2 |AckOH lAdlOlj . and the 
particle spectra [KH03 well. These models predict that vi will rise monotoni- 
cally with p?. It is expected however, that for particles with large pt the model 
assumptions will break down. Measurements using charged particles and our 
measurements using Kg and A + A do indeed indicate a saturation of vi — well 
below the hydrodynamic calculations — at intermediate p?- The nuclear modifica- 
tion factor Rep for charged particles also shows a large pr-independent suppres- 
sion of particle production at p? > 4.5 GeV/c. Models based on parton energy 
loss |GVW01I IShu02| and transport opacity |MG02| have been discussed in rela- 
tion to the saturation and centrality dependence of v-i at intermediate and high 
Pt- The authors of Reference pTVW02j propose that the saturation of charged 
particle v 2 could be a consequence of the transition from soft to hard production 
processes occurring at different p^-scales for pions and protons. 

We show here that the extent of the pt region where hydrodynamic like pro- 
cesses (or other soft processes) dominate the spectrum is particle species depen- 
dent. In Chapter EJ we reported the measurement of t>2 and Rcp for Kg and 
A + A. These measurements show that for px up to 3 GeV/c A + A v 2 continues 
to rise — similar to hydrodynamic model calculations — while the K s v 2 saturates 
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at pt ~ 1.5 GeV/c. We also find that there is a particle-type dependence to the 
onset of the suppression — as measured by Rcp — of K s and A + A production. In 
addition, for each particle the lower px bound of the suppressed region of Rcp 
coincides with the lower bound of the saturated region of vi- 

It has been suggested that if a partonic state exists prior to hadronization, the 
process of particle formation at intermediate p?, by string fragmentation, par- 
ton fragmentation |LK02b| or quark coalescence |LK02al lEMUHl IH Y031 IHKL03al 
FMN03b, FMN03a , may lead to a dependence of v% and Raa on particle type. 
In this case, it is possible that these measurements will provide information on 
the existence and nature of an early partonic state. In this chapter we investigate 
the interplay between the apparently soft and hard components of the Kg and 
A + A spectrum and explore possible sources for the particle-type dependence 
of our measurements. We begin with a brief description of hydrodynamic and 
energy loss models and we compare their predictions to our data. In Section I5~31 
we use hydrodynamical inspired and pQCD inspired fits to estimate the pt where 
soft, hydrodynamical-type processes are no longer appreciable. 

5.1 Describing Heavy-Ion Collisions with Hydrodynamics 

A system can be described within a hydrodynamical formalism when the time 
scales of its microscopic processes are sufficiently smaller than the time scale 
for its macroscopic evolution. For heavy-ion collisions this means the time be- 
tween interactions amongst the constituents — partonic and/or hadronic — must 
be much smaller than the lifetime of the system. When this condition is met, 
the constituents can interact enough times to equilibrate. The space-time evo- 
lution of the system can then be described in the framework of relativistic fluid 
dynamics. The equations of motion are derived from the conservation of energy 
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and momentum d^T^ = 0. The energy-momentum tensor T^ u in the ideal fluid 
approximation is given by: 



= (e + pK 




(5.1) 



where e, p, and are respectively the energy density, pressure, and four velocity. 

Prior to the onset of local equilibrium, the hydrodynamical equations are 
invalid. As such, they can only describe heavy-ion collisions from an initial time 
To until the time when interaction rates in the system become too small and local 
thermal equilibrium can no longer be maintained (i.e the freeze-out time). The 
initialization of the hydrodynamic evolution requires that the pre-thermalization 
stage be modelled so that the initial conditions can be estimated. Given a set 
of initial conditions and the hydrodynamic equations, all that remains to be 
specified is the nuclear equation-of-state EOS which relates the thermodynamic 
quantities of the system. The EOS can be modelled or calculated using lattice 
QCD. Finding the EOS that governs nuclear matter at high temperature and 
density is the primary objective of heavy-ion physics and hydrodynamical model 
calculations may provide insight into its form. 

The hydrodynamic equations for an ideal fluid describe the velocity and pres- 
sure fields for thermalized fluid elements. The Cooper-Frye formula |CF74j is 
used to calculate the momentum distribution for the hadrons created from the 
fluid elements on the freeze-out hyper-surface S: 



where d{ is a degeneracy factor, ^ are the chemical potentials for the hadrons, 
p^ are their four momentum, and do^ is the outward normal vector on S. Hy- 



drii di f P^da^ 



(5.2) 



d 3 p (2vr) 3 A exp[(p^ - m)/T*\ =F 1 ' 
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drodynamic models have been particularly successful in reproducing the mass 
dependence of v 2 at low pr- Figure 15.11 shows t>2 versus pr (left) and the inte- 
grated t>2 versus particle mass for identified particles (right) compared to hydro- 
dynamical model calculations. In these models the increase in the integrated v 2 
with mass is a consequence of the collective motion of the fluid elements built up 
as the system interacts and expands. When particles with more mass freeze-out 
from fluid elements flowing with a given velocity, they will carry greater mo- 
menta. In this way, an anisotropic collective flow velocity at freeze-out leads to 
an increase of the integrated v 2 with particle mass. This increase and the mass 
dependence of the differential v 2 both indicate that a significant collective motion 
is established — perhaps early in the collision — and that the source eccentricity is 
efficiently transferred to momentum-space anisotropy. 




Figure 5.1: Left: Differential v% at low px for pions, kaons, protons, and lamb- 
das at mid-rapidity from Au+Au collisions at y/s NN = 200 GeV. The charged pion, 
charged kaon, and the proton i>2 were measured by the PHENIX collaboration |Esu03j . 
Right: Integrated 1)2 versus mass compared to hydrodynamic calculations with different 
freeze-out temperatures. The integrated V2 is calculated by weighting the measured V2 
with the particle yield extracted from fits to the spectra. 



The results of a full hydrodynamical model calculation can be approximated 
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with a simple analytical model: the blast wave model HKHOlj. The equations 



of the blast wave model describe particle emission from a thin cylindrical shell of 
thermalized matter with temperature T. The particle production is approximated 
by a boosted Boltzmann distribution so that the particle spectra can be calculated 
from the equations: 



d 3 



n 



dp 



3 



2tt 

rf0 s iT 1 (/5(0 s ))e a( ^ )cos( ^-^ ) , (5.3) 



a(x,p) 



Tsinh(p(a;)) ' 
rriT 

Tcosh(p(x)) ' 



where (f) p = tan -1 (p y /p x ), and <p s = t&n~ 1 (y/x) are respectively the momentum 
and coordinate space azimuthal freeze-out angles, K\ is the modified Bessel's 
function and p(x) is the transverse flow rapidity. In Section 1573*1 we'll use the blast 
wave equations to parameterize the soft part of the K s and A + A spectrum 1 . 

5.2 Energy Loss 

For p+p collisions, when scatterings involve sufficiently large momentum transfer, 
pQCD calculations describe hadron production well (see Figure lTTiT) . To relate the 
partonic and hadronic invariant cross-sections {E a ^^- and E^^- respectively) 
it is assumed the calculation of the partonic cross-section and the hadronization 
process for production of hadron h with momentum p can be factorized |Owe87j : 



1 The author thanks F. Retiere for his help with these fits. 
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d 3 n h ^ f 1 dz d 3 n a 



dp 3 ^Jo z 2 dp A a 
The probability that parton a with momentum p a fragments into hadron h with 
momentum p = z x p a is expressed in terms of fragmentation functions D a ^h(z)- 
The fragmentation functions are typically taken to be universal: Once measured 
they can be used to describe hadron production for other hard processes. 

For heavy-ion collisions however, neither the validity of factorization nor 
the universality of the fragmentation functions can be assumed a priori. Fast 
partons — presumably produced from hard interactions between two colliding 
nuclei — may need to traverse hot, dense matter before escaping from the sys- 
tem. Over two decades ago Bjorken estimated that these secondarily produced 
quarks and gluons could lose tens of GeV of their initial transverse momentum 
via elastic scattering with quanta in the medium |Bjo82| . In this publication 
Bjorken also proposed what would later be called surface emission as a signature 
of jet-quenching. 

An interesting signature may be events in which the hard collision 
occurs near the edge of the overlap region, with one jet escaping 
without absorption and the other fully absorbed. 

Recent observations at RHIC confirm that in central Au+Au collisions, while 
near- angle jet-like correlations exist, away-side jet-like correlations are suppressed 
Adl 03f] IAda03bj . Jet-quenching has been studied extensively throughout the 



last decade (see for example References [TG9H IWG921 IGW941 IZak971 |BDM97, 



IGLVOOl IWieOUl IBDM01I IMul03| ) and although energy loss from elastic scattering 
has been shown to be small, radiative energy loss may be large |GP90l IBD P95 , 
WG92J. The magnitude of the energy loss is thought to depend on the gluon 



density of the medium being traversed, and should therefore be sensitive to the 
creation of hot, dense, and perhaps deconfined matter. 

Measurements of the suppression of neutral pions and charged hadrons in cen- 
tral Au+Au collisions confirm many of the expectations of energy loss jAdcOll 
Adl 02bt IAda03bt IAdl03f| , with the suppression accompanied by the disappear- 
ance of back-to-back jet-like correlations [Adl03a, Adl03bJ. It's possible (even 
likely) however, that the interactions of the fast partons with the matter, will 
induce not only energy loss but also changes to the hadronization process (DH77, 
IRHM79[lAltMlAMH01[imM02l^ 

RS03, FMN03a]. In this case, even if factorization still proves to be valid, we 
can't assume the fragmentation functions measured in e + + er collisions will be 
relevant to Au+Au collisions. Measurements of the production of identified par- 
ticles are needed to study the possible evolution of the fragmentation functions 
with system size and to understand the processes that may govern those changes. 
In so doing, it may be possible to not only reach a better understanding of heavy- 
ion collisions but to also develop a much deeper understanding of hadronization 
in general. This would constitute a major advance relevant to all particle physics, 
where previously hadronization has only been dealt with phenomenologically. 

Surface emission has been discussed in relation to the large, pr-independent v 2 
measured for charged hadrons |Shu02j and the large pr-independent suppression 
of charged particle production at high p?. Within this scenario, more energy loss 
will lead to a larger v-i and a greater suppression of particle production. This is 
inconsistent, however, with our measurements of v<i and Rcp for kaons and A + A 
at intermediate px'- we find that kaons have a smaller v 2 but a larger suppression. 
These calculations do not, however, account for how the process of hadronization 
may change the observed v-i. 
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Figure 5.2: Comparison of Kg and A + A with expectations from surface emission 
models (see text). 

Figure O shows v 2 for K Q S and A + A from the 30-70%, 5-30%, and 0- 
5% centrality intervals along with calculations of v 2 from two surface emission 
scenarios: "W.S. matched" uses a Woods-Saxon distributions in the calculation 
of the nuclear overlap function and matches the energy loss to the observed 
suppression of neutral pion production [Jia02l while "H.S. maximum" uses an 
analytic function representing pure surface emission — infinite energy loss — from a 
hard sphere overlap geometry lVol03j . The models are not in good agreement with 
our measurements. We do not, however, rule out the surface emission scenario 
for pt above 4.5 GeV/c. We also note that for pt from 1.5-4.5 GeV/c, it may be 
that kaons are produced predominantly from hard processes in a surface volume 
while ( ant i-) lamb das in the same range are produced by soft processes. The 
magnitude of the kaon v 2 in this region, however, is much larger than would be 
expected from the surface emission model with a realistic nuclear overlap density. 
Stronger conclusions on this point can be drawn from a larger data sample and 
more extensive studies of possible non-flow contributions to identified particle v 2 . 
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In the case that particles are produced from hard processes — either near the 
surface or throughout the entire volume — their spectra are expected to be well 
represented by a power-law function [Alb901 . In Section l5~!?l we use the power-law 
function 

£?~CU + ^r (5-5) 
dp 6 p 

to parameterize the hard part of the Kg and A + A spectrum. 
5.3 Transverse Momentum Regimes 



In what follows, we study the two component nature of the Kg and A + A 
spectra. We seek to delineate the boundary between the soft and hard region of 
the spectra. 

In Reference jSor03j we stressed the effectiveness of the combination of Rep 
and Vi for mapping out the transition between the region dominated by soft pro- 
cesses and the region dominated by hard processes. More detailed calculations 
making use of these ideas in Reference |HN03j are in relatively good agreement 
with the available data. We use hydrodynamical model inspired blast wave func- 
tions and pQCD motivated power-law functions to fit the Kg and A + A spectra 
with the intention of extracting the value of pt where the soft-to-hard crossover 
occurs. Figure 15.31 shows blast wave and power-law fits to the Kg and A + A 
spectra for three centrality intervals. The fit parameters are listed in Table l5~Tl 

Radial flow, as established in hydrodynamical models for example, can lead 
to a mass dependence in the soft-to-hard px crossover PT,cross- Reference (HN03 
uses the pt where the yields from hard and soft production processes become 
equivalent to define Pt, cross- Using hydrodynamic model calculations for the soft 
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Transverse Momentum p,- (GeV/c) 



Figure 5.3: Blast wave and power-law fits to the Kg and A + A spectra. Fitting 
parameters are listed in Table I3TTI 

contribution, and a pQCD parton model — incorporating energy loss, gluon shad- 
owing, and initial state rescattering — for the hard component, they calculate 
PT,cross = 1-8, 2.7, and 3.7 GeV/c for pions, kaons, and protons respectively. 
From our blast- wave and power- law fits, however, judging by the applicability 
of the fit functions we find that pr,cross for kaons in central events is closer to 
1.8 GeV/c, and consistent with the pion PT,cross apparent in the pion spectrum 



Fit 

Parameters 


0-5% 40-60% 60-80% 


A(A) 

0-5% 40-60% 60-80% 


C 

Po (GeV/c) 
a 


330.11 
3.280 
16.43 


3.6785 
2.786 
13.96 


0.0829 
2.717 
13.41 


40.914 
2.12e+5 
5.58e+5 


0.9832 
10.38 
34.25 


0.0424 
5.691 
21.91 


T (MeV) 
Po 


140.7 
0.686 


179.1 
0.499 


175.5 
0.489 


136.3 
0.866 


258.7 
0.499 


223.2 
0.412 



Table 5.1: Fitting parameters for power-law (first three rows) and blast wave (last 
two rows) parameterizations of the Kg and A + A spectra. 



92 



measured by PHENIX |Ad.l03dj . Given the pt reach of our measurements it's 
difficult to determine pr,cross for lambdas from spectra fits alone. We can, how- 
ever, conclude that the lambda PT,cross is greater than 3 GeV/c: within a range 
that would be consistent with the proton pT,cr OSS quoted in Reference [HN03 . 
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Figure 5.4: The Rcp of K° and A + A with blast wave curves extracted from fits to 
the particles spectra. 



The pt where the A + A Rcp (0-5%/60-80%) begins to decrease in Figure 
(bottom panel) provides clearer evidence for the cross-over and suggest that for 
( ant i-) lamb das PT,cross ~ 3.5 GeV/c. We note, however, that the blast wave fits 
to the A + A spectra for the 0-5% and 40-60% centrality intervals are in good 
agreement throughout the fit range so that we cannot exclude the possibility that 
the blast wave parameterization describes lambda production well throughout the 
measured pr- This is also apparent in the blast wave curves for Rcp using those 
same centralities (the top panel of Figure I5~4"j) . To determine whether or not 
this indicates that the assumptions of the hydrodynamical models are valid from 
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0-60% centrality but break down in the 80-100% interval will require further 
study. 

The similarity of the Pt,ctoss f° r pions and kaons may indicate that the number 
of quarks rather than the mass of the hadron provides the relevant scale for the 
transition between regions of predominantly soft and predominantly hard pro- 
duction processes. This empirical observation is consistent with a picture where 
at intermediate px, hadronization occurs through the coalescence of co-moving 
partons. The recent observation that Rcp for the meson (m=1.019 GeV/c 2 ) is 
closer to the pion and kaon Rcp than the proton or lambda Rcp MY03 supports 
this picture. 
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Figure 5.5: The V2 and Rcp of Kg and A + A. The vertical lines emphasize the 
correlation between the saturation of V2 and the decline of Rcp- 



To emphasize the correlation between the behavior of identified particle v% 
and identified particle Rcp we plot Vi and Rcp together in Figure 1531 Although 
Rcp depends only on the yield in the central and peripheral bins, and the v% 
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in Figure 15.51 is from a minimum-bias centrality interval, the two parameters 
may still be intimately related. Differential elliptic flow v^ipr) measures the 
difference between the pt spectrum of particles emitted in the direction of the 
reaction plane (in-plane) to that of particles emitted perpendicular the reaction 
plane (out-of-plane). In hydrodynamic models, we expect the pressure gradient 
to be larger in the in-plane direction than the out-of-plane direction. Since the 
average pressure gradient in central collisions is expected to be larger than in 
peripheral collisions in a the hydrodynamical picture the ratio Rqp and v 2 should 
have a similar px and particle-type dependencies. Coalescence or recombination 
models should also give a similar correlation between u 2 and Rcp- Both measures 
reflect relative probabilities for forming hadrons: The probabilities depend on the 
phase-space density of partons, and the phase-space density varies with centrality 
and azimuthal angle. In Appendix we introduce a formalism more suited to 
studying the centrality dependence of particle yields in heavy-ion collisions across 
a broad pt range and for a variety of system sizes. 

5.4 Initial State Effects 

Nuclear enhancement of hadron yields (relative to number-of-binary-collisions 
scaling) for lower beam energies has been observed at intermediate pt in p+A 
collisions — with a larger enhancement for baryons than mesons [Str92j. The 
empirical observation of an enhancement in the scaled yield is known as the 
Cronin effect |(Jro73llAnt79j . The Cronin effect is generally attributed to multiple 
scatterings between the projectile partons and the cold nuclear matter in the 
target |Acc02j . The ratio of the scaled yields in d+Au and p+p collisions RdAu 
has been measured for neutral pions |Adl03c| and inclusive charged hadrons at 
RHIC [Ada 03a| IBac03at IAdl03ct IArs03| . At intermediate px, the neutral pion 
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RdAu is consistent with one (no Cronin effect) or a small enhancement. The 
inclusive charged hadron RdAu shows an enhancement of roughly 35% however, 
indicating the presence of a particle-type dependent Cronin effect at RHIC energy. 
The existence of an enhancement has not been established for kaons, lambdas, 
or antilambdas, but the Cronin effect cannot be ruled out as the origin of the 
dependence of Rep on particle type. 

Theoretical calculations, such as those involving initial parton scatterings off 
cold nuclear matter (e.g. |LP83j ). don't reproduce the particle- type dependence 
of the enhancement factor observed in p+A collisions. The inadequate particle 
dependence in these calculations may arise from the fact that these models only 
deal with initial parton scatterings while the observed hadrons are formed at the 
late stage of the collision. It is argued in Reference fAG03| . that the fragmenta- 
tion process can distort the features of the parton level Cronin effect. As such, the 
strong particle-type dependence in R p (d)A may indicate a nuclear modification of 
the parton fragmentation into baryons and mesons or alternatively, the presence 
in p(d)+A collisions of a multi-parton particle formation mechanism such as co- 
alescence [GKL03aJ or recombination [FMN03aJ. These mechanisms are beyond 
the framework of many existing theoretical models for the Cronin effect. 



5.5 Hadronization of Dense Matter 



The varied initial conditions in heavy-ion collisions may provide important 
phenomenological clues to the nature of hadron formation — particularly baryon 
formation. Observing how the presence of dense nuclear matter influences hadroniza- 
tion will almost certainly help clarify by what processes three quarks can come 
together to form a baryon. 
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Modification of the fragmentation functions D a -+h(z) has been offered as a 
strategy to account for the influence of the surrounding matter on the fragmen- 
tation process [GWOOl IWGOlj . When this strategy leads to a rescaling of z it 
should affect all hadrons in a similar way |FMN03aj — a behavior that is incon- 
sistent with the measurements presented in this thesis. To account for species 
dependence observed in t> 2 and Rcp a new understanding of hadronization may 
be necessary. In addition, we note that in the complex and highly interacting 
systems created in heavy-ion collisions the assumption of factorization implicit 
in Equation 15.41 may no longer be valid. 

Although the observables of heavy-ion collisions — Rcp, particle ratios, and V2 
for example — and their variation with particle-type challenge current theoretical 
models, they also hint at possible resolutions. 

We note for example, that the absence of a significant suppression with re- 
spect to binary scaling of the A + A yield at intermediate pr in central Au+Au 
collisions may indicate the presence of dynamics beyond parton energy loss and 
standard fragmentation. The larger A + A Rcp at intermediate pt means that 
the (anti-)lambda yield increases with the parton density of the collision fireball 
much faster than the meson yield. The rate of increase of the proton yield |Adl03ej 
and the multi-strange baryon yield (S + 5) |Lon03| are found to be similar to 
that of (anti-)lambda. In addition, at intermediate pr, the centrality depen- 
dence of the 0-meson production is more similar to that of kaons than (anti- 
lambdas |MY03j . Stronger dependence on parton density for baryon production 
is naturally expected from multi-parton production mechanisms such as gluon 
junctions |V(t99| . quark coalescence |MV03j . or recombination |FMN03bJ. 

As discussed in Section 15. 3 1 it is apparently the number of constituent quarks 
in a hadron that predominantly determines the characteristics of it's spectra at 
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intermediate pt- Figure ESI shows v-i of Kg and A + A as a function of pt, where 
the V2 and px values have been scaled by the number of constituent quarks (n). 
While i>2 is significantly different for K s and A + A, within errors, 1*2/11 vs Pt/ 11 is 
the same for both species above Pt/ 11 ~ 0.8 GeV/c. This behavior is consistent 
with a scenario where hadrons at intermediate pt are formed from bulk partonic 
matter by coalescence of co-moving quarks: in this case ^/n vspr/n reveals the 
momentum space azimuthal anisotropy that partons develop from the collision 
ellipsoid, e.g. Reference |MV03| . Hadronization by coalescence and the large 
partonic anisotropy subsequently inferred from our empirical observations would 
both argue strongly for the existence of a strongly interacting early partonic stage. 



Figure 5.6: The V2 parameter for Kg and A + A scaled by the number of constituent 
quarks (n) and plotted versus pr/ n - 

The suggestive scaling behavior if Figure lo~o1 can be tested by including mea- 
surements of t>2 for other identified particles that extend into the intermediate 
Pt region. Figure IK77I shows t>2/n versus pr/ri for all identified particles currently 



available from RHIC experiments |Ksu03| [UB03J. At low p T (p T < 1 GeV/c), 
where hydrodynamic calculations were already seen to reproduce the measured 
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t>2 well, and at px > 6 GeV/c, the signatures of coalescence are not expected to 
be prominent. The deviation of the scaled pion v-i may prove to be problematic 
or it may just reflect the break-down of coalescence at low px- Otherwise, the 
results in Figure 15.71 are consistent with the expectations of constituent-quark- 
number scaling and provide strong evidence for coalescence and the existence of 
a quark-gluon plasma in the early stage of the collision system. 
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Figure 5.7: The V2 parameter for K® and A + A scaled by the number of constituent 
quarks (n) and plotted versus pr/n. 



For nxpT,hard, where pr,hard is a momentum for which the spectra of the 
underlying partons follows a power-law, the yield of hadrons from coalescence 
of n quarks will no longer dominate the yield of hadrons from fragmentation of 
partons with momentum pT,hard/z. Based on this picture, we would expect the 
V2 of all hadrons at high px {pr > 6) to take on the unsealed value of the parton 
v 2 . Assuming the W2/ 11 values from Figure 15.71 reflect the parton v 2 we would 
conclude that at high px (where we expect v^ eson = v ^ ar y° n — v ^ uark ^ V2 Q f all 
light-flavored particles (particles having only u, d, or s constituent quarks) will 
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take on the same value V2,hard ~ 0.07. This assumes a saturated parton v% at high 
Pt'- an assumption that is consistent with expectations from transport [MG02J, 
and surface emission [Shu02j models. 

5.6 Conclusions 

In summary, we have reported the measurement of v 2 and Rqp up to p? of 
6.0 GeV/c for kaons and A + A from Au+Au collisions at y/s NN = 200 GeV. At 
low Pt, hydrodynamic model calculations agree well with i> 2 f° r K s anc ^ A + A. 
At intermediate p?, however, hydrodynamics no longer describes the particle pro- 
duction well. For K s , v 2 saturates earlier and at a lower value than for A + A. 
In addition, Rcp shows that the kaon yield in central collisions is suppressed 
more than the (anti-)lambda yield. At intermediate pt, the A + A yield in cen- 
tral Au+Au collisions is close to expectations from binary scaling of peripheral 
Au+Au collisions. At high p?, the Rcp of K s , A + A and charged hadrons are ap- 
proaching the same value. The measured features in the kaon and (anti-)lambda 
t>2 and Rcp may indicate the presence of dynamics beyond the framework of 
parton energy loss followed by fragmentation. The particle- and pr-dependence 
of v 2 and Rcp, particularly at intermediate pt, provides a unique means to in- 
vestigate the anisotropy and hadronization of the bulk dense matter formed in 
nucleus- nucleus collisions at RHIC. 

In Au+Au collisions that copiously populate phase space, it would be naive 
to expect a single parton description of hadronization to remain valid. Our 
measurements verify that multi-parton dynamics significantly alter production 
mechanisms: giving rise to a number-of-constituent-quark dependence for particle 
spectra, Rcp, and v 2 at intermediate pj>. From the mass dependence of v 2 at low 
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Pt and the partonic v-i inferred from the scaled v 2 at intermediate px it appears 
that the existence of a thermalized partonic state at RHIC is not just likely but 
perhaps unavoidable. 



5.7 Future Directions 



The future direction of our studies at RHIC are clear: verification of the cre- 
ation of a QGP, then characterization of the QGP. The measurement of identified 
particle v 2 and Rcp will play a prominent role in both of these endeavors. The full 
implication of the measurements presented in this thesis will become clearer when 
t> 2 and R C p have been measured well into the hard region (p T ~ 8 GeV/c) for 
7T°, Kg, and <fi mesons and for p, A, S, and Q baryons. Several of these measure- 
ments are already available and most of the others will become available in the 
near future. The measurements already made argue strongly for the existence of 
a thermalized partonic state that at least partly hadronizes by coalescence. With 
the larger data sets collected during future RHIC runs v 2 should be measured for 
identified particles using a cumulant analysis. The cumulant analysis can easily 
be adapted for statistically identified particles when only one particle candidate 
in a given p T interval is considered from each event. 

Once a suite of identified meson and baryon measurements has been made to 
sufficiently high pt, the first priority should be to conduct a system-size scan. 
RHIC is already operating at its top energy for Au+Au collisions and running 
with lower energy will drastically reduce high pr particle yields — essentially elim- 
inating several of our most promising QGP signatures. Instead, it will be more 
fruitful to study collisions with fewer participants. This will provide three excit- 
ing opportunities: A chance to search for a region where the QGP state is turned 
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on or off, a chance to run at even higher energy with the existing RHIC facility 
(Si+Si top energy would be ^/s NN = 250 GeV) and the chance to search for 
rare/exotic states (i.e. glue-balls, penta-quarks, etc.) in collisions that generate 
less combinatorial background. The reduction in the combinatorial background 
will benefit many identified particle measurements. Cu+Cu collisions may be 
preferable since they should still allow reliable measurements of v 2 while probing 
smaller system sizes. 
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APPENDIX A 



Collision Geometry and the Source Eccentricity 



Here we provide information about the geometry and coordinate systems of rel- 
ativistic heavy-ion collisions. 




Figure A.l: Coordinate system for the transverse plane in heavy-ion collisions. Con- 
tours of the overlap density are also shown. Here <fi is the azimuthal angle with respect 
to the y-axis. Also shown are generalized coordinates el and e~2 used to calculate the 
fourth-harmonic eccentricity where the (y-)e~2-axis is the only uniquely defined direction 
in the transverse plane. 
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Figure IA~T1 gives the coordinates in the azimuthal plane of the collisions. The 
beam direction (z-axis) is in or out of the page as indicated by the ® or © 
respectively. The d/-)e~2-axis is uniquely defined (using the right-hand-rule) by 
the collision axis and the directions of the colliding beams. When calculating 
t>2 = (cos(20)), <p is defined with respect to the reaction plane (the x-axis). For 
our eccentricity calculations we will define (f) with respect to the y-axis. In this 
way, v 2 and the 2 nd -harmonic eccentricity will have the same sign. 

The density of Au nucleus is represented by a Woods-Saxon distribution, 

1 . . 

PA = PO— , r -R A s ( AA ) 

where r = y/s 2 + z 2 , R A = 1.12 x A 1/3 , with p = 0.159 GeV/fm 3 and f = 0.535 
fm. The thickness function of nucleus A (T^) is the nuclear density integrated 
over the z direction: T A (s) = f dzpA(a,s)- The density of nucleons participating 
in the collision (wounded nucleons) is given by: 

n part {x, y) = T A x (1 - e ~ TBr7 ™) +T B x(l- e" T ^), (A.2) 

where for collisions at y/s NN = 200 GeV we use a^N ~ 42. The density of binary 
nucleon-nucleon collisions is taken as: nu n (x,y) = a^^ x T4 x Tg. In the plots 
that follow, we use n part as the density in the overlap region. 

To better understand the initial spatial asymmetry of the collision regions in 
Figure lA~2l we show calculations of eccentricities defined for even harmonics (the 
odd harmonic eccentricities are all zero). The most commonly used definition of 
the second harmonic eccentricity is €2 = (y 2 — x 2 )/(y 2 + x 2 ). In the figure we 
also include calculations of eccentricities from 2(cos(n0)) and a calculation of the 
4* ft -harmonic eccentricity using the appropriate coordinate system as shown in 
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Figure A. 2: Left: Analytic calculations of the even ni/i-harmonic eccentricities ver- 
sus impact parameter. Right: Relative magnitudes of harmonics for various impact 
parameters. For more central collisions the higher harmonics fall off more quickly. 

Figure IA.1I We note that the magnitude of higher harmonic eccentricities falls 
off more quickly in central collisions than peripheral collisions (panel b). For 
peripheral collisions the higher harmonics are appreciably larger. The different 
methods for calculating the eccentricity do not yield the same result. In Fig- 
ure IA.3I we show the eccentricity from the more standard calculation. We also 
show the ratio of the 2 nd - and 4 i/l -harmonic eccentricities versus impact parameter 
(panel b). 
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Figure A. 3: Left: Analytic calculations of the 2 nd -harmonic and 4 th -harmonic eccen- 
tricities. Right: The ratio of the second and fourth harmonic eccentricities. 
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In Figure IA.4I (a) we plot profiles (along the x and y axis) of the overlap 
density for the impact parameters listed in panel b. In Figure lA~4l (b) we plot the 
gradient of the profiles in panel a. In a hydrodynamical picture — once a system 
has thermalized — the density changes will generate pressure gradients that drive 
the flow of matter and induce collective motion. We note that Figure IA.4I (b) 
suggests that the greatest initial pressure gradients will be found in the in-plane 
direction of collisions having impact parameter 6 = 6 fm. This impact parameter 
corresponds roughly to collisions within the 10-20% centrality interval. 




Figure A. 4: Left: Collision overlap densities calculated for various impact parameters 
(listed in panel b) and using a Woods-Saxon nuclear density profile and a wounded 
nucleon model. Right: The gradients of the density profiles in panel a. We take the 
magnitude of the gradient to estimate of the initial pressures in the system. 
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In Figure IAT5I we show the geometries corresponding to the centrality intervals 
used in our v 2 measurements. The mean eccentricities and impact parameters in 
this figure were calculated as in Reference |Adl02cj 1 . 
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0-5% 



5-30% 
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Figure A. 5: Overlap densities, mean eccentricities, and mean impact parameters for 
the centrality intervals used for our v<i analysis. 



The author thanks A. Tang for his help with these calculations. 
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Figure A. 6: Configurations for the collision centrality intervals used to calculate 
Rcp- 

In Figure IA.6I we show the geometry of the centrality intervals used in the 
calculation of Rcp- 
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APPENDIX B 



Kinematic Variables 

The azimuthal components of a particles momentum p x and p y are used to define 
its transverse momentum, 

Pt = \jvl + Pi- (B.l) 
The transverse mass/energy of a particle having mass m is 

m T = \Jp 2 T + m % ( B - 2 ) 

so that the transverse kinetic energy of the particle is tut — tuq. The transverse 
kinetic energy is commonly used in place of px- 

In the lab frame, the azimuthal angle of a particles momentum is simply 
(fiiab = tan _1 (pj,/pa;). When the reaction plane can be measured the azimuthal 
angle can be measured with respect to the transverse angle of the reaction plane 

<P = <Plal> - ^R P . (B.3) 

This is the shown in Figure IA.1I where the reaction plane is aligned with the 
x-axis. 

With the transverse coordinates defined all that remains is to define a longi- 
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tudinal variable. The rapidity y is defined as 



y 




E + p z 
E-Vz 



) 



(B.4) 



and is boost invariant. In the case that the momentum of a particle is known but 
not its energy — typically because it's mass is unknown — then the pseudo-rapidity 
rj = — In tan (cos~ 1 (p z /p)/2) can be used. For p >> m, rj ~ y. 

The differential cross-section for particle production, is found by counting the 
number of particles d 3 n produced in a phase space element. It's advantageous to 
define the phase-space element to be Lorentz invariant. The usual choice is the 
element dp x dp y dp z /E, so that the invariant cross-section is 

d^n d 3 n _ , 

V = ^Wv' (R5) 

where we've expressed the cross-section in terms of the variables defined previ- 
ously. 
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APPENDIX C 



An Improved Formalism for Studying the 
System Size Dependence of Nucleus-Nucleus 

Collisions 

Up to now the centrality dependence of particle yields in heavy-ion collision has 
primarily been studied either by plotting the spectra from different centralities 
in the same panel or by forming the ratio 

R , s [dn[ (N b inaryCfor)p ntral (Q ^ 

CPKPT> [rfn/(N binary dp r )] pcriphcra1 ' 

This ratio is formed in order to test how different high p T particle production 
in central collisions is from peripheral collisions where yields are expected to 
scale by the number of binary collisions. Here we suggest another formalism that 
we believe is better suited to the study of the centrality dependence of particle 
production in a variety of heavy-ion collision systems and across the entire pr 
range. 

Particle yields scale predominantly with the number of participating nucleons 
Np art , not Nbinary and Np art is less model dependent than Nonary F° r this rea- 
son we choose to consider the yield of a particle Y from collisions with impact 
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parameter b scaled by N part : 

y(b) = y(6)/N part . (C.2) 
We can write y(b) in terms of a Taylor expansion around a given impact parameter 

y{b) = y<M (' + ml {b - w + ^k> - + ■ ■ ■) • (a3) 

We label the terms si = ^§g, s 2 = , etc. so that s n = (n _ 1) 1 !y(6o) gg. In 

this case if the yield scales with N part then s n = for all ra. The term si represents 
the linear dependence of yip) and s 2 represents the quadratic dependence. 

These measures have several benefits over the scaled ratio Rcp- First we note 
that if one of the centrality parameters s n is dominant then coalescence scaling 
would give s™ eson (p T ) = 2 x s q n aark (p T /2) and s b n aryon (p T ) = 3 x s q n aark (p T /3). Even 
if several parameters must be taken into account trivial scaling rules relating 
meson and baryon s parameters can still be attained. Second, we can easily 
define a minimum-bias s n measurement or an s n measurement within a given 
centrality so that the centrality dependence at different intervals and for different 
systems can be more easily compared. For example, we can expand yip) about 
bo = b max /2.0 and calculate the minimum-bias si = 12 ^~ b ») ( w here the scaled 
yield y has been used to calculate the mean). These variables may be particularly 
valuable for studying changes to the systematic variation of the yield with system 
size. One would hope to see a notable change in the s parameters as the system 
size increases to a size necessary to form a QGP. 

A correction to the s n variables can be introduced to account for the impact 
parameter resolution. 
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APPENDIX D 



The STAR Collaboration 



J. Adams 3 , C. Adler 11 , Z. Ahammed 24 , C. Allgower 12 , J. Amonett 14 , 
B.D. Anderson 14 , M. Anderson 5 , D. Arkhipkin 10 , G.S. Averichev 9 , J. Balewski 12 , 

O. Barannikova 9 ' 24 , L.S. Barnby 14 , J. Baudot 13 , S. Bekele 21 , V.V. Belaga 9 , 
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